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I.  Introduction 


This  is  the  final  report  on  work  performed  under  the  support  of  AFOSR 
Grant  F49620-85-K-0021  during  the  period  1  Oct.  1885  to  30  Sept.  1887.  The 
investigation  covered  the  generation  of  radiation  from  relativistic  electron 
beams.  The  scope  and  findings  of  the  investigation  are  summarized  in  section 
II.  Further  details  of  this  work,  in  reviewed  papers,  preprints  and  reports, 
are  appended. 


II.  Summary  of  Work  Accomplished 

During  the  funding  period  1  Oct.  198S  to  30  Sept.  1987,  we  carried  out 
AFOSR  supported  research  into  mode  competiton  in  the  Lowbitron  and  efficiency 
enhancement  through  magnetic  field  profiling;  a  novel  method  of  gain  and 
efficiency  enhancement  in  the  CARM  (Electron  Cyclotron  Autoresonance  Maser) 
which  makes  use  of  a  transverse  helical  wiggler  to  inhibit  parasitic 
autoresonant  acceleration  of  electrons  by  the  electromagnetic  wave;  absolute 
instablility  and  tunability  of  the  CARM  through  varying  the  solenoidal  field; 
and  cyclotron  effects  in  an  FEL  (Free  Electron  Laser)  with  an  axial  guide 
magnetic  field.  We  also  investigated  the  technical  specifications  for  a  CARM 
oscillator  in  the  megawatt  power  range,  with  the  aim  of  designing  a  cost 
effective  experiment  utilising  an  existing  Magnetron  In^^ticn  Gun  and  a 
gyrotron  magnet,  which  will  serve  to  verify  the  basic  physics  of  the  CARM  and 
demonstrate  its  viability.  Extensive  use  of  computer  simulations  using 
relativistic,  electromagnetic  PIC  codes  was  made  in  these  studies. 


II. 1  Mode  Competition  In  the  Lowbltron 


The  lowbitron  is  a  Doppler-shift  dominated  cyclotron  maser  with  a 
longitudinal  wiggler  field  such  that  the  net  magnetic  field  is  of  the  form.  ? 
=z(Bc+BySlnkwz) ,  where  By  and  Bc  are  respectively,  the  wiggler  and  solenoldal 
fields.1’2  The  lowbitron  has  certain  attractive  features.  It  is  capable  of 
generating  coherent  radiation  at  a  frequency  w=(k+kw)vz+Qc  «  2r2 ( kwvz+Qc ) , 
where  vz  is  the  beam  velocity  along  the  solenoldal  field  and  £lc  is  the 
relativistic  cyclotron  frequency.  This  output  frequency  is  higher  than  that 
produced  by  a  free  electron  laser  with  the  same  wiggler  modulation,  kw,  or 
that  produced  by  a  Doppler-shift  dominated  cyclotron  maser  with  the  same  Qc. 
Furthermore,  it  has  the  practical  advantage  in  that  the  longitudinal  wiggler 
is  easier  to  produce,  and  the  wiggler  modulation  easier  to  change,  than  is 
the  case  with  the  transverse  wiggler  in  an  FEL.  However,  it  had  been  thought 
that  the  efficiency  of  this  device  would  be  low. 

A  l2/2~d  particle  code**  has  been  used  to  study  the  lowbitron.  The 
simulations  (Appendix  1)  show  that  the  cyclotron  mode  w=kvz+ftc,  and  the 
lowbitron  mode  u=(k+kw)vz+fic  can  be  simultaneously  excited  by  the  electron 
beam.  Competition  between  these  two  modes  lead  to  a  finite  threshold  in 
Bw/Bc  for  efficient  operation.  Figure  1  shows  the  efficiency  and  growth  rate 
as  a  function  of  Bw.  The  growth  rate  Initially  decreases  as  Bw  increases, 
until  a  threshold  at  Bw~Bc/2  is  exceeded,  after  which  the  growth  rate 
Increases  linearly  with  Bw  and  the  lowbitron  mode  dominates.  Efficiency  as 
high  as  252  is  attained  in  the  simulations. 


The  efficiency  of  the  lowbitron  can  be  further  increased  with 
appropriate  tapering  of  the  solenoidal  field  to  prolong  the  deceleration 
phase  of  the  electron-wave  interaction.  The  dynamic  phase  shift  of  the 


lowbitron  is  given  by  A$  =  [  <k+k  >Av  -  (fl  /y>Ay  It,  where  t  is  the 
electron  transit  time  in  the  interaction  region.  As  the  rate  of  change  of 
phase  due  to  v and  Y  are  not,  in  general,  equal,  the  electron  eventually 
drifts  into  the  acceleration  phase  with  respect  to  the  wave,  and  so  ceases 
the  transfer  of  energy  to  the  wave.  To  prolong  the  deceleration  phase,  the 
solenoidal  field  can  be  tapered  so  that  A$  remains  small.  We  find  that  it 
is  possible  to  raise  the  untapered  field  efficiency  of  25 7.  to  more  than  40 7. 
by  increasing  the  solenoidal  field  with  time. 


I I . 3  Enhancement  of  Gain  and  Efficiency  in  a  CARM 

The  Doppler-shift  dominated  cyclotron  autoresonance  maser  (CARM)  is  an 
attractive  source  of  coherent  radiation  in  the  millimeter  and  submillimeter 
wavelength  range. ’  It  possesses  a  favourable  ouput  frequency  scaling  with 
electron  beam  energy,  as  well  as  the  potential  for  high  efficiency.  The  wave 
frequency  for  resonant  wave-electron  interaction  is  w=kvz+  ftc,  which  scales 
approximately  with  the  beam  kinetic  energy  as  w*y2qCj  y,  being  the  Lorentz 
factor  of  the  beam.  For  relativistic  phase  velocity  Vp=w/k*c,  the  change  in 
Qc  is  largely  compensated  by  the  change  in  the  Doppler  term  when  an  electron 
interacts  with  the  wave.  This  phase-locking  mechanism  between  electron  and 
wave  provides  the  basis  for  high  efficiency.  However,  one  drawback  is  that 
the  gain  decreases  as  Vp-»c.  Recalling  that  an  electron  can  lose  as  well  as 


gain  energy  In  Its  Interaction  with  the  wave,  the  decrease  in  gain  can  be 
viewed  as  increasing  competition  between  the  processes  of  stimulated  electron 
emission  and  absorption. 


We  have  found  that  both  the  gain  and  efficiency  can  be  significantly 
enhanced  by  adding  a  transverse  helical  wiggler  field  to  the  CARM  in  the 
interaction  region.  The  wiggler  field  is  that  produced  by  a  bifilar  current 
winding  used  in  present  free  electron  lasers,  given  by  the  model  equation 
?w=B„<xCosk„z  +  ySinkwz).  <The  net  static  magnetic  field  configuration  is 
the  same  as  that  in  an  FEL  with  a  guide  magnetic  field,  but  unlike  the  FEL, 
the  electron  beam  is  injected  into  the  wiggler  region  with  a  finite 
transverse  velocity. )  The  wiggler  is  used  to  effect  selective  inhibition  of 
autoresonance  acceleration  of  electrons  (or  stimulated  absorption)  by  the 
wave.  In  the  combined  wiggler  and  uniform  axial  magnetic  field,  the 
unperturbed  electron  velocity  has  an  oscillating  component,  vz»v0+pCos(Qc- 
ky v o ) ,  where  vQ  is  the  mean  value  and  p« (ftc“kwvo)_1  ls  small,  awav  from 
gyroresonance  <Qc=kwv0).  By  choosing  appropriate  values  for  Bw  and  kw,  it 
may  be  arranged  such  that  electrons  that  undergo  autoresonance  acceleration 
will  approach  gyrosonance,  acquiring  large  oscillations  in  vz  which  quicklv 
detune  the  electrons  from  resonance  with  the  wave,  and  thereby,  terminating 
the  acceleration  process.  On  the  other  hand,  electrons  that  lose  energy  to 
the  wave  will  suffer  less  perturbation  in  vz,  allowing  them  to  stay  locked  in 
phase  with  the  wave.  Further  details  of  this  enhancement  mechanism  are  given 
in  Appendix  3.  One  and  2/g-d  particle  simulations  show  that  significant 
enhancement  in  gain  and  efficiency  might  be  achieved.  Figure  2  shows  an 
example  of  the  efficiency  and  gain  in  a  CARM  with  (solid  dots),  and  without 
(hollow  dots)  a  wiggler  field,  as  a  function  of  the  wave  phase  velocity 
|3p=w/kc . 


I I . 4  Stability  and  Tunabllitv  of  a  CARM  Amplifier 

Frequency  tunablllty  of  a  CARM  amplifier  may  be  achieved  by  varying 
either  the  electron  beam  energy  or  the  magnetic  field.  The  viability  of  the 
latter  scheme  has  been  Investigated  using  theoretical  analysis  and  computer 
simulations.  Stability  consideration  dictates  an  upper  limit  to  the  magnetic 
field  that  can  be  used  to  upshift  the  frequency  in  a  CARM.  If  the  cyclotron 
frequency  is  too  high,  the  cyclotron  mode  can  excite  a  low  frequency  back 
propagating  wave,  as  well  as  the  desired  high  frequency  forward  propagating 
wave,  resulting  in  an  absolute  instability  (oscillation).  One  way  of 
suppressing  the  oscillation  is  to  use  a  resistivelv  coated  waveguide.  An 
example,  using  a  TEoi  mode,  and  an  electron  beam  with  the  parameters  (767keV, 
2.SkA,  vi/c=0.4)  is  considered.  Analytic  theory  developed  by  Briggs^  and  Chu 
et  al6  is  used  to  predict  the  threshold  for  the  absolute  instability.  It  is 
found  that  a  40-60%  variation  of  the  magnetic  field  can  be  achieved  for  a 
maximum  wall  resistivity,  parameterized  by  6=0.01,  where  6  is  the  ratio  of 
wall  skin  depth  to  wall  radius.  A  40%  variation  in  the  magnetic  field, 
however,  can  lead  to  a  frequency  tuning  range  of  over  100%.  Computer 
simulations  for  these  beam  parameters  show  linear  growth  rates  in  good 
agreement  with  theory.  Efficiency  in  excess  of  20%  can  be  achived  before 
wave  saturation  by  electron  phase  trapping.  The  effects  of  beam  momentum 
spread  are  also  studied  in  the  computer  simulations.  They  indicate  that  high 
beam  quality  is  crucial  at  high  operating  frequency.  The  CARM  ceases  to  be 
an  attractive  device  for  beam  momentum  spread  greater  than  3%.  Further 
details  are  given  in  Appendix  4. 
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I I . S  CARM  Oscillator  Design  and  Analysis 

Investigations  into  the  CARM  to  date  have  been  confined  mostly  to 
theoretical  studies.  In  light  of  the  favorable  theoretical  expectations,  it 
is  of  pressing  interest  to  design  a  specific  experiment  which  would  allow 
verification  of  our  understanding  of  the  basic  physics  of  the  CARM  and 
demonstrate  its  viability  as  a  practical  device.  With  these  alms  in  mind,  we 
have  initiated  a  study  into  the  technical  specifications  for  a  cost  effective 
CARM  oscillator.  The  CARM  oscillator  would  operate  in  the  megawatt  power 
range  and  would  make  use  of  an  existi  ig  Magnetron  Gun  and  gyrotron  magnet. 

The  study  analysed  design  of  the  Bragg  reflector  and  resonant  cavity.  Mode 
competition  with  the  backward  wave  is  expected  to  be  a  important  problem  and 
this  was  analysed  for  specific  oscillator  configurations.  Based  on  the 
analysis,  a  tentative  design  for  a  CARM  oscillator  is  as  follows: 


voltage .  400  kV 

current . 25  A 

v  / v . 0.92 

i  z 


operating  frequency  .... 

. 150 

GHz 

cavity  0  . 

. 376 

cavity  length  . 

. 5.1 

cm 

cavity  radius  . 

. 29 

cm 

axial  magnetic  field  .  .  . 

. 37 

kG 

beam  fill  factor  . 

.  407. 

I I . 6  Hybrid  Cyclotron  Instability  In  an  FEL  with  Axial  Guide  Field 

From  simulations  using  a  l^^-d  particle  code,  it  is  found  that  a  hybrid 
cyclotron  instability  can  be  excited  at  ^(k+ZkyJvQ-Oc ,  with  relatively  large 
growth  rate  and  high  efficiency.  This  instability  may  be  explained  by  the 
following  simple  picture.  In  the  combined  wlggler  and  axial  magnetic  field, 
the  beam  electrons  have  a  natural  frequency  of  axial  oscillations,  given 
approximately  by  w'=lflc-  kwv0  1  (asuming  a  right  helical  wiggler).  An 
electromagnetic  wave  with  frequency  w=(k+2kw)v0  will  produce  an  axial 
ponderomotive  force  at  the  resonant  frequency  of  electron  axial  oscillations. 
The  resonant  ponderomotive  force  can  lead  to  large  amplitude  oscillations  in 
the  electron  axial  velocities,  which,  in  turn,  is  coupled  through  the  wiggler 
field  to  give  rise  to  a  transverse  current  which  reinforces  the  wave  at 
u»=<k+2kw)v0-  flc.  The  existence  of  a  natural  frequency  of  oscillation  in  the 
electron  motion  was  first  pointed  out  in  an  earlier  analytic  work  by 
Frledland  and  Hlrshfield.7  They  were  primarily  concerned  with  using  this 
effect  to  enhance  the  free  electron  laser  gain  very  close  to  gyroresonance , 
or  the  unstable  beam  propagation  condition  (at  gyresonance,  the  frequencies 
of  tl  hybrid  and  FEL  mode  (u)=(k+kw)v0>  become  the  same).  Subsequent 
authors®'^  have  also  ignored  the  existence  of  a  separate  hybrid  mode. 

Results  from  our  simulations  (see  Appendix  6)  show  that  the  hybrid  mode  can 
be  distinguish  from  the  FEL  mode.  Also,  and  more  importantly,  an 
autoresonance  effect  appear  to  occur  under  certain  circumstances,  giving  rise 
to  high  efficiency  <>30Z)  for  the  hybrid  mode.  This  hybrid  mode  may  lead  to 
new  versatility  for  the  free  electron  laser. 
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Competition  of  electron-cyclotron  maser  and  free-electron  laser  modes 
with  combined  solenoidal  and  longitudinal  wiggler  fields 

A.  T.  Lin  and  Chih-Chien  Lin 

Department  of Physics,  University  of  California  at  Los  Angeles,  Los  Angeles,  California  90024 
(Received  26  December  1985;  accepted  7  February  1986) 

A  relativistic  electron  beam  with  a  finite  transverse  dc  momentum  (Bt  =  l/y0)  passing  through  a 
region  of  combined  uniform  solenoidal  and  longitudinal  wiggler  magnetic  fields  is  observed  to 
convert  25%  of  its  kinetic  energy  into  coherent  radiation  at  frequency  to  =  yl  (kw  VQ  +  tle/y0)  if 
the  phase  velocity  of  the  generated  wave  is  slightly  above  the  speed  of  light  In  this  situation,  the 
bunching*  of  the  slow  electron-cyclotron  mode  and  free-electron  laser  modes  with  combined 
solenoidal  and  longitudinal  wiggler  fields  (lowbitron)  are  observed  to  compensate  each  other, 
which  gives  rise  to  a  finite  threshold  for  lowbitron  operation.  In  order  to  attain  high  efficiency,  the 
wiggler  strength  of  a  lowbitron  must  substantially  exceed  the  threshold. 


In  the  millimeter  and  submillimeter  wavelength  ranges, 
the  free-electron  devices  being  actively  pursued  at  many  in¬ 
stitutions  include  the  gyrotron,1  Doppler-shift-dominated 
cyclotron  maser,1"1  and  free-electron  laser.4"1  In  all  these 
devices,  the  energy  exchange  process  is  one  of  converting 
relativistic  electron  beam  kinetic  energy  into  coherent  radi¬ 
ation  energy  by  slowing  the  electrons  down  through  electro¬ 
magnetic  forces.  In  order  for  electrons  to  give  up  their  ener¬ 
gy  efficiently,  they  must  be  bunched  in  the  decelerating 
phase  of  the  generated  wave  and  remain  in  this  favorable 
phase  as  long  as  possible.  The  bunching  force  can  come  ei¬ 
ther  from  a  fast  wave  or  a  slow  wave.  In  these  wavelength 
ranges,  most  experimental  efforts  are  concentrated  on  gyro 
trons  and  free-electron  lasers.  However,  a  very  attractive 
alternative  to  these  devices  is  the  Doppler-shift-dominated 
cyclotron  maser.  This  mechanism  for  generating  microwave 
radiation  possesses  the  potential  of  achieving  high  efficiency, 
as  well  as  a  favorable  output  frequency  scaling  with  the  elec¬ 
tron  beam  energy.  Recently  it  has  been  shown  theoretical¬ 
ly6,7  that  a  relativistic  electron  beam  with  a  finite  transverse 
dc  momentum  passing  through  a  region  of  combined  uni¬ 
form  solenoidal  and  longitudinal  wiggler  magnetic  fields 
(lowbitron)  is  capable  of  generating  coherent  radiation  at 
frequency  aa2yi(kw  V0  +  Clc/y0),  where  fle  and  V0  are, 
respectively,  the  solenoidal  field  electron-cyclotron  frequen¬ 
cy  and  electron  axial  velocity,  kw  is  the  wiggler  wavenumber, 
and  y0  is  the  electron  relativistic  factor.  The  output  frequen¬ 
cy  is  tunable  by  varying  the  electron  beam  energy  and  is 
higher  than  the  frequency  from  either  a  free-electron  laser 
with  the  same  wiggler  wave  number  or  a  Doppler-shift- 
dominated  cyclotron  maser  with  the  same  solenoidal  field. 
However,  the  energy  conversion  mechanism  is  thought  to  be 
a  low-efficiency  process.  On  the  other  hand,  it  has  been  dem¬ 
onstrated1  through  computer  simulations  that  a  Doppler- 
shift-dominated  cyclotron  maser  operated  in  an  autoreson¬ 
ant  region  can  convert  25%  of  the  beam  kinetic  energy  into 
coherent  radiation.  In  this  paper,  we  carry  out  computer 
simulations  to  show  that  using  the  magnetic  configuration 
described  above  also  has  the  potential  to  achieve  high  effi¬ 
ciency. 

Consider  a  thin  pencil  beam  of  relativistic  electrons  that 
traverses  on  axis  through  a  magnetic  field  region  of  the  form 


B  =  i(B0  +  Bu,  sin*„z).  (1) 

The  strength  of  the  solenoidal  (B0)  and  wiggler  (BK )  mag¬ 
netic  fields  can  be  comparable  in  experiment.  The  electron 
transverse  motion  under  the  influence  of  the  magnetic  field 
expressed  by  Eq.  ( 1 )  is  a  superposition  of  Larmor  and 
wiggler  motion  so  that  a  transverse  electromagnetic  wave 
(eo,k)  becomes  resonant  with  the  electron  motion  if  the  fol¬ 
lowing  condition  is  satisfied: 

to  =  (*r  +  k„  )v t  +  Clc/y.  (2) 

In  emitting  a  transverse  electromagnetic  wave  propagating 
along  the  axial  direction,  an  electron  resulting  from  the  re¬ 
coil  of  wave  momentum  will  lose  both  its  energy  and  axial 
momentum  such  that 


where  Z  =  m^fy  and  is  the  phase  velocity  of  the  wave. 
Replacing  P,  by  m0yv ,  in  Eq.  (3)  yields 

Ly^yJ^W'/d-B+Bo),  (4) 

where  3  —  v/c.  In  order  for  an  electron  interacting  with  a 
wave  to  lose  its  energy  as  well  as  its  axial  velocity,  K0  must  be 
less  than  c2/^ .  In  this  situation,  the  phase  shifts  caused  by 
the  Doppler  and  cyclotron  terms  will  compensate  each  oth¬ 
er.  In  doing  so,  the  electron  can  remain  in  resonance  with  the 
wave  for  a  long  time,  and  it  converts  most  of  its  kinetic  ener¬ 
gy  into  radiation.  The  compensation  becomes  exact  if  the 
phase  velocity  is  close  to  the  speed  of  light  and  the  resonance 
condition  is  fulfilled  initially.  This  can  be  seen  by  examining 
the  expression  for  the  dynamic  phase  shift, 

=  -  [(*,+*„)*». +n£A(l/y)]r,  (5) 

where  T  is  the  electron  transit  time  in  the  interaction  space. 
Substituting  Eqs.  (2)  and  (4)  into  Eq.  (5)  yields 


Since  the  second  term  of  Eq.  (6)  must  be  negative  for  an 
electron  to  lose  both  its  energy  and  axial  velocity  at  the  same 
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time,  the  phase  velocity  of  the  generated  wave  must  be 
greater  than  the  speed  of  light  for  the  dynamic  phase  shift  to 
remain  small.  In  a  lowbitron,  the  condition  for  to  be  zero 
is  no  longer  =  c  but  is  Fj*  >c.  This  is  to  compensate  the 
phase  shift  arising  from  an  electron  executing  wiggler  mo¬ 
tion  in  a  longitudinal  wiggler  field. 

Without  the  longitudinal  wiggler  field,  the  relativistic 
electron  beam  with  a  finite  dc  transverse  momentum  propa¬ 
gating  through  a  region  of  uniform  solenoidal  field  is  already 
unstable  with  respect  to  electron-cyclotron  maser  instabil¬ 
ity.  The  introduction  of  a  longitudinal  wiggler  field  will  be 
considered  as  a  perturbation  to  the  original  unstable  system. 
The  dispersion  relation  of  the  electron-cyclotron  mode  can 
be  expressed  as* 

^■^4— - 

Yo \v-  k,v0  -  *Vy0 

l(a  -  k,v o  -  Clc/Yo)2J 

where  cof  is  the  plasma  frequency  of  the  electron  beam. 
There  are  two  types  of  instabilities  with  different  bunching 
mechanisms.  In  the  regime  of  Fpb>c  (fast  wave),  the 
bunching  mechanism  depends  on  azimuthal  electron  bunch¬ 
ing  caused  by  the  relativistic  change  in  electron  mass  when 
electrons  interact  with  the  transverse  electric  field  of  the 
wave  (the  so-called  negative  mass  effect).  Since  an  electron 
gyrates  faster  in  a  magnetic  field  as  it  loses  energy,  the  elec¬ 
trons  initially  tend  to  bunch  in  phases  where  the  electric 
force  ahead  is  accelerating  and  that  behind  is  decelerating.  If 
the  Doppler-shifted  wave  frequency  is  slightly  larger  than 
the  relativistic  electron-cyclotron  frequency,  the  electron 
bunches  slip  behind  the  wave  in  phase  and  on  average  trans¬ 
fer  their  energy  to  the  wave.  On  the  other  hand,  for  the  case 
of  K,*  <c  (slow  wave),  the  longitudinal  Lorentz  force  pro¬ 
duced  by  the  cross  product  of  the  electron  dc  transverse  ve¬ 
locity  and  the  magnetic  field  of  the  wave  provides  the  neces¬ 
sary  bunching.  The  electrons  initially  tend  to  bunch  in 
regions  where  the  force  ahead  is  decelerating  and  the  force 
behind  is  accelerating.  If  the  Doppler-shifted  wave  frequen¬ 
cy  is  slightly  smaller  than  the  relativistic  electron-cyclotron 
frequency,  the  electron  bunches  gain  on  the  wave  in  phase 
and  on  average  transfer  their  energy  to  the  electromagnetic 
wave.  Therefore,  the  electron  bunching  for  the  fast-wave  and 
the  slow-wave  branches  should  take  place  180*  out  of  phase 
in  the  electron  phase  space.  According  to  Eq.  ( 2 ) ,  the  impo¬ 
sition  of  a  longitudinal  wiggler  field  will  shift  the  phase  of 
electron  bunching  by  roughly  Fe¬ 
in  order  to  investigate  the  effects  of  a  longitudinal 
wiggler  on  the  gain  and  efficiency  of  an  electron-cyclotron 
maser,  a  one-and-two-halves-dimensional  relativistic  elec¬ 
tromagnetic  particle  code9  is  used.  The  following  parameters 
are  chosen:  y0  =  2- 5,  PL  =  l/yo,  <op  =  0.35<uc,  kwVu 
=  0.81a>c,  and  <uc  =  fle/y0.  In  order  to  attain  high  effi¬ 
ciency,  the  electrons  must  remain  in  resonance  with  the 
wave  for  a  long  time,  so  a  single  mode  operation  is  essential. 
This  can  be  accomplished  by  using  either  a  Bragg’s  resonator 
for  an  oscillator  configuration  or  injecting  a  signal  with  in¬ 
tensity  well  above  the  noise  level  in  an  amplifier  configure- 
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tion.  Only  the  results  of  single  k  mode  simulation  will  be 
given.  The  phase  velocity  of  the  generated  wave  and  the 
wiggler  strength  are  varied  to  determine  the  parameters  for 
achieving  high  efficiency.  The  electron  phase-space  struc¬ 
ture  [(y-  1)  vs  ^  =  o>rt -  k,zt  -  tan ~'(.P,/P,)]  at  the 
time  when  substantial  electron  bunching  has  occurred  is  dis¬ 
played  in  Fig.  1.  The  results  of  Bw  =  0  are  shown  in  Fig.  1(a) 
for  a  fast  wave  (V^  =  1.037c)  and  Fig.  1(b)  for  a  slow  wave 
( Fp,,  —  0.93c).  It  is  very  clear  that  the  phase  of  electron 
bunching  for  these  two  branches  differs  by  180*. 

By  increasing  the  amplitude  of  the  wiggler  field,  the  low¬ 
bitron  mode  within,  =  (k  +  kw )  V0  +  Clc/y0  will  gradually 
become  dominant  over  the  electron-cyclotron  maser  mode 
with  Wo,,  =  kV0  +  Clc/Yo-  In  order  to  investigate  the  effects 
of  their  competition  on  the  gain  and  efficiency,  a  series  of 
simulations  varying  the  wiggler  field  from  Bw  =  0  to  Bw 
—  1.2B0  has  been  caried  out.  The  wavenumber  is  chosen  to 
be  kc  =  [Fig.  1(b)]  so  that  FJ*  sza^/k  =  1.004c, 
which  meets  the  condition  of  autoresonance  of  Eq.  (6).  The 
time  evolution  of  the  generated  wave  energy  for  two  typical 
cases  (Bw  =  0  and  Bw  =  B0)  is  shown  in  Fig.  2(a).  For  the 
case  of  Bw  =  0,  the  wave  energy  exhibits  an  exponential 
growth  at  early  time  with  a  growth  rate  of  to,  —  0.05  lcoc, 
which  is  close  to  estimation  from  Eq.  (7)  (a>t  =  0.05 5oc ) 
and  saturates  at  a  level  that  is  less  than  1%  of  the  beam 
kinetic  energy.  The  wave  energy  of  the  case  of  Bw  —  B0  in¬ 
creases  with  a  growth  rate  of  co,  —  0.036a>c ,  but  saturates  at  a 
much  higher  level,  which  is  about  25%  of  the  beam  kinetic 
energy.  The  growth  rate  and  efficiency  as  a  function  of  the 
wiggler  strength  are  given,  respectively,  in  Figs.  2(b)  and 
2(c).  The  simulation  results  reveal  that  the  growth  rate  first 


T  i  5 


FIG  1.  The  phase  space  of  electron  -cyclotron  maser  modes  at  v,t  -  100; 
(a)  fast-wave  branch  ( ^  —  1.037c)  and  (b)  slow-wave  branch 
(K„  -0.93c). 
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FIG.  2.  The  effects  of  longitudinal  wiggler  on  the  electron-cyclotron  maser 
instability;  (a)  the  time  evolution  of  the  generated  wave  energy,  (b)  growth 
rate  versus  wiggler  strength,  and  (c)  efficiency  as  a  function  of  wiggler 
strength. 


decreases  as  the  wiggler  strength  is  increased  until  the 
wiggler  strength  exceeds  a  threshold,  which  in  this  case  is 
roughly  0. 5 B0.  For  Bw  >  0.6.00,  the  growth  rate  increases  lin¬ 
early  with  the  wiggler  strength  and  a  lowbitron  begins  to  be 
able  to  attain  high  efficiency. 

In  order  to  understand  the  threshold  formation  of  a  low¬ 
bitron,  consider  the  wiggler  field  as  a  perturbation  to  the 
electron-cyclotron  maser  mode.  The  increase  in  wiggler  field 
gradually  changes  the  unstable  frequency  from  the  cyclotron 
maser  mode  tc  the  lowbitron  mode  ain, .  This  is  demon¬ 
strated  in  Fig.  3,  which  shows  the  frequency  spectrum  as  a 
function  of  wiggler  strength.  At  Bw  =  0,  the  frequency  spec¬ 
trum  peaks  at  the  cyclotron  maser  frequency 
co  —  Oca  =  9  04<u,.  With  Bw  =*  0.5 B&  two  discrete  peaks, 
one  at  a>  -  <wOT  and  the  other  one  at  a  =  -  9.84<ut ,  ap¬ 

pear  in  the  frequency  spectrum.  When  Bm  is  increased  to  Bv 
the  peak  at  co  —  ty,*  becomes  dominant.  According  to  Eq. 
(2),  there  are  two  discrete  frequencies  (eo^  and  con )  for 
each  k  mode.  Similarly,  in  an  amplifier  system,  an  input  coQ 
will  excite  two  modes  with  discrete  wavenumbers.  These  two 
modes  with  a  frequency  difference  of  kw  V0  tend  to  interfere 
with  each  other,  which  gives  rise  to  a  finite  threshold  of  the 
wiggler  field  for  a  lowbitron.  The  time  history  of  electron 
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FIG.  3.  The  frequency  spectrum  u  s 
function  of  wiggler  strength;  (a) 
B.-O,  (b)ff.«  0.5*0,  (O  B, 
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FIG.  4.  The  time  history  of  phase  space  for  a  lowbitron  with  *„ 
<u, t  -  0,  (b)  <u,i  -  100,  (c)  a>,» «•  120. 
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phase  space  for  a  lowbitron  with  Bm=B0  is  shown  in  Fig.  4. 
The  phase  of  the  maximum  electron  bunching  differs  from 
Fig.  1(b)  by  an  amount  of  roughly  kw  V0  [Fig.  4(b)],  and 
Fig.  4(c)  indicates  that  phase  trapping  is  the  dominant  satu¬ 
ration  mechanism. 

In  conclusion,  the  competition  between  the  electron-cy¬ 
clotron  maser  and  lowbitron  modes  is  believed  to  give  rise  to 
a  finite  threshold  for  a  lowbitron.  By  choosing  the  param¬ 
eters  such  that  Fj*  ><-,  a  lowbitron  is  capable  of  attaining 
high  efficiency  of  the  order  of  25%. 
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A  HIGHLY  EFFICIENT  ELECTRON  CYCLOTRON  MASER  IMCRSEO  IN  A  LONGITUDINAL  HIGGLER 
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Abstract 

A  relativistic  electron  beam  with  a  finite 
transverse  dc  aoaentun  (flj.  "  1/t0)  passing  through 
s  region  of  coablned  uniform  solenoldal  and  longi¬ 
tudinal  wlggler  magnetic  field  la  observed  to 
convert  25Z  of  Its  kinetic  energy  Into  coherent 

2 

radiation  at  frequency  w  •  y  (k  V  -*■  Q/v) 

o  wo  co 

If  the  phase  velocity  of  the  generated  wave  la 
slightly  above  the  speed  of  light.  In  order  to 
attain  high  efficiency,  the  wlggler  strength  of  a 
lowbltron  must  be  comparable  to  the  solenoldal 
field  strength.  Furthermore,  the  efficiency  can 
be  Increased  to  about  42Z  If  the  solenoldal  field 
is  properly  profiled. 

Introduction 

In  the  millimeter  and  submllllmeter  wave¬ 
length  ranges,  the  free  electron  devices  being 
actively  pursued  at  many  Institutions  Include 
gyrotron  (1],  Doppler-shifted  dominated  cyclotron 
maser  (2]  (3],  and  free  electron  laser  (4]  [3],  In 
all  these  devices,  the  energy  exchange  process  Is 
one  of  converting  relativistic  electron  bean 
kinetic  energy  Into  coherent  radiation  energy  by 
slowing  the  electrons  down  through  electromagnetic 
forces.  In  order  for  electrons  to  give  up  their 
*h*rgy  efficiently,  they  must  be  bunched  In  the 
decelerating  phaae  of  the  generated  wave  and 
remain  in  this  favorable  phase  as  long  as  possi¬ 
ble,  The  bunching  force  can  come  either  from  a 
fast  wave  or  a  slow  wave.  In  these  wavelength 
ranges,  most  experlasntal  efforts  are  concentrated 
on  gyrotroes  and  free  electron  lasers.  However,  a 
*•*7  attractive  alternative  to  these  devices  Is 
the  Doppler-sblft  dominated  cyclotron  maser.  This 
.mechanism  of  generating  microwave  radiation  pos- 
*••••0  the  potential  of  achieving  high  efficiency, 
*•  well  as  a  favorable  output  frequency  scaling 
with  the  electron  beaa  energy.  Recently  It  has 
keen  shown  theoretically  (6]  17]  that  a  relativls- 
electron  beaa  with  a  finite  transverse  dc  mom¬ 
entum  passing  through  a  region  of  combined  uniform 
adenoidal  and  longitudinal  wlggler  magnetic 
fields  (lowbltron)  Is  capable  of  generating  co- 


2  a. 

herent  radiation  at  frequency  w  *  2y  (k-V.  ♦  — ) 

o  "  To 

where  Qe  and  V0  are  respectively  the  solenoldal 
field  electron  cyclotron  frequency  and  electron 
axial  velocity,  k^  is  the  wlggler  wave  number,  and 
Y0  Is  the  electron  relativistic  factor.  The  output 
frequency  la  tuneable  by  varying  the  electron  beaa 
energy  and  Is  higher  than  the  frequency  from 
either  a  free  electron  laser  with  the  same  wlggler 
wave  number  or  a  Doppler-shlft  dominated  cyclotron 
maser  with  the  same  solenoldal  field.  However, 
the  energy  conversion  mechanism  Is  thought  to  be  a 
low  efficient  process.  On  the  other  hand,  It  has 
been  demonstrated  [3]  through  computer  simulations 
that  a  Doppler-shlft  dominated  cyclotron  maser 
operated  In  an  auto-resonant  region  can  convert 
23Z  of  the  beam  kinetic  energy  Into  coherent 
radiation.  In  this  paper,  we  carry  out  computer 
simulations  to  show  that  using  the  magnetic  con¬ 
figuration  described  above  also  has  the  potential 
to  achieve  very  high  efficiency  If  the  solenoldal 
magnetic  field  Is  properly  profiled. 

Simla  Electron  Calculation 

Consider  s  thin  pencil  beaa  of  relativistic 
electrons  traverses  on  axis  through  a  magnetic 
field  region  of  the  form 

I  -  x  (Bo  +  Bw  sin  kw*)  (1) 

Ths  strength  of  the  solenoldal  (B0)  and  wlggler 
(By)  magnetic  field  can  be  comparable  In  real 
experiment.  The  electron  transverse  motion  under 
the  influence  of  the  magnetic  field  expressed  by 
Bq.  (1)  is  s  superposition  of  Larmor  and  wlggler 
motion  so  that  a  transverse  electromagnetic  wave 
(w,  k)  becomes  resonant  with  the  electron  motion 
If  ths  following  condition  Is  satisfied: 

m  ■  fk  ♦  k  )  v  ♦  — §  (2) 

'  x  w'  x  Y 

In  emitting  s  transverse  electromagnetic  wave 
propagating  along  axial  direction  an  electron  due 


to  recoil  of  wave 


tntum  will  lose  both  Its 


energy  end  axial  momentum  such  that 


J  2  2 

WZ  -  k£C  « 


(7) 


<3) 

where  E  -  ■0c2T  and  vph  !•  the  phase  velocity  of 
th«  wav*.  Replacing  P,  by  In  8q  •  (3)  yialda 

T  6  .48 

«-iVH  •  <‘> 

ph  o 

♦ 

what*  I  -  -i.  In  order  for  an  alactroo  Interacting 
c 

with  a  wava  to  loaa  lta  energy  aa  well^aa  lta 

axial  valoclty,  Va  auat  ba  laaa  than  | — .  In  thla 

vph 

situation,  tha  phasa  ahlft  caused  by  tha  Dopplar 
and  cyclotron  tarma  will  coapanaata  aach  othar. 

In  doing  so,  tha  alactroo  can  raaaln  In  raaonanc* 
with  tha  wava  for  a  long  tla*  and  convart*  aost  of 
lta  kinetic  anargy  Into  radiation.  Tha  co a pant  a  - 
tlon  becomes  exact  If  tha  phaaa  valoclty  la  eloaa 
to  tha  spaad  of  light  and  tha  raaonanc*  condition 
la  fulfilled  Initially.  This  can  ba  saan  by  exau- 
log  tha  exprasslon  for  the  dynaalc  phasa  ahlft 

48  -  -((k£  ♦  kw)  4v£  ♦  Qe4  (I)]  T  ,  <5) 

where  T  la  tha  alactroo  transit  tlaa  In  the  Inter¬ 
action  space.  Substituting  Eqs.  (2)  and  (4)  Into 
Eq.  (3)  yallda 

f?c  t , _ 1  ■>  _  *w«  r,  _  ,  a 


“-*1 


T.  (6) 


Since  tha  second  tens  of  Eq.  (6)  aust  ba  negative 
for  an  electron  to  lose  both  Its  anargy  and  axial 
valoclty  at  tha  a aaa  tlaa,  tha  phaaa  valoclty  of 
tha  generated  wava  auat  ba  greater  than  tha  spaed 
of  light  for  tha  dynaalc  phaaa  shift  to  raaaln 
saall.  In  a  lowbltron,  tha  condition  for  48  to  ba 
aero  Is  no  longer  “  c  but  is  Vph  >  c.  This  Is 
to  coapanaata  tha  phasa  shift  arising  froa  an 
electron  executing  wlgglar  notion  In  a  longitudi¬ 
nal  wtggler  field. 

Electron  Cyclotron  Haaar 

Without  tha  longitudinal  wlgglar  field,  tha 
relativistic  electron  beaa  with  a  finite  dc 
transverse  aoaantua  propagating  through  a  region 
of  unifora  solanoldal  field  la  already  unstable 
with  raapact  to  electron  cyclotron  aasar 
Instability.  Tha  Introduction  of  a  longitudinal 
wlgglar  field  will  ba  considered  as  a  perturbation 
to  the  original  unstable  aysteu.  The  dispersion 
relation  of  electron  cyclotron  node  can  ba 
axprassad  aa  (8] 


w  -  k  v 


2  2  2  2 
8.  (<*  -  k  c  ) 


p  f  ~  Vo  “1  7  , 

To  [u-kv  '  7(«-kv  f1 

*  o  Vo  *  o  Y0  ‘ 

where  *p  la  the  plasaa  frequency  of  the  electron 
beaa.  There  are  two  types  of  Instabilities  with 
different  bunching  aechanlsa.  In  tha  raglae  of 
Vph  >  c  (fast  wave),  tha  bunching  aachanlsa 
depends  on  aalauthal  electron  bunching  caused  by 
the  ralatlvlatlc  change  In  electron  aaaa  whan 
electrons  Interact  with  the  transverse  electric 
field  of  the  wava  (tha  so-called  negative  aass 
effect).  Since  aa  electron  gyrates  faster  in  a 
aagnatlc  field  aa  It  loses  energy,  the  electrons 
Initially  tend  to  bunch  In  phases  where  tha 
electric  force  ahead  Is  accelerating  and  that 
behind  la  decelerating.  If  the  Doppler-shifted 
wave  frequency  la  slightly  larger  than  the  relati¬ 
vistic  electron  cyclotron  frequency,  the  electron 
bunches  slip  behind  the  wava  In  phase  and  on 
average  transfer  their  energy  to  the  wave.  On  the 
other  hand,  for  the  case  of  Vpj,  <  c  (slow  wave), 
the  longitudinal  Lorenta  force  produced  by  tha 
crosa  product  of  tha  electron  dc  transverse  valoc¬ 
lty  and  tha  nagnetlc  field  of  the  wave  provides 
the  necessary  bunching.  The  electrons  Initially 
tends  to  bunch  In  regions  where  the  force  ahead  Is 
daceleratlng  and  the  force  behind  Is  accelerating. 
If  the  Doppler-shifted  wave  frequency  la  slightly 
saaller  than  tha  relativistic  electron  cyclotron 
frequency,  the  electron  bunches  gain  on  the  wave 
In  phasa  and  on  average  transfer  their  energy  to 
the  electromagnetic  wave.  Therefore,  the  electron 
bunching  for  the  fast-wave  and  tha  slow-wave  bran¬ 
ches  should  take  place  180*  out  of  phase  In  the 
electron  phase  space.  According  to  Eq.  (2),  tha 
laposltlon  of  a  longitudinal  wlgglar  field  will 
ehlft  the  phase  of  electron  bunching  by  roughly 

kwV 

In  order  to  Investigate  the  effects  of  a 
longitudinal  wlgglar  on  the  gain  and  efficiency  of 
an  electron  cyclotron  naser,  an  on*  and  two  halves 
dlaanslonal  relativistic  elect roaagnetlc  particle 
cod*  [9]  Is  used.  The  following  paraaetere  are 


chosen:  Y  -  2.5,  8.  -  i-, 

*  -  *  To 


0.35  w  ,  k  V 
c  wo 


e 

0.81  u  and  u  «  — .  In  ordar  to  attain  high 
c  C  T<| 

efficiency,  the  electrons  aust  remain  In  resonance 
with  Che  wave  for  a  long  tin*  so  a  single  nod* 


operation  la  aaaantlal.  Thin  can  bo  accomplished 
by  ualnf  either  a  Bragg 'a  reaouator  for  an  oacll- 
lator  configuration  or  injecting  a  signal  with 
intensity  well  above  the  nolae  level  in  an  amplif¬ 
ier  configuration.  Only  the  result*  of  single  k 
mode  simulation  will  be  given.  The  phase  velocity 
of  the  generated  wave  and  the  wlggler  strength  are 
varied  to  determine  the  parameters  for  achieving 
high  efficiency.  The  electron  phase  space  atruct- 

p 

ore  ((T  -  1)  versus  4  -  w_t  -  k-ij  -  tan-1  at 

rx 

the  time  when  substantial  electron  bunching  has 
occured,  as  well  as  the  bunching  force,  are  dis¬ 
played  la  Pig.  1.  The  results  of  iy  ■  0  are  shown 
In  Pig.  lb  for  fast  wave  (Vpj,  ■  1.037c)  and  Pig. 
lc  for  slow  wave  (Vph  ■  0.93c).  It  la  very  clear 
that  the  phase  of  electron  bunching  for  these  two 
branches  differs  by  180  degrees. 

It  has  been  shown  In  (10)  that  with  By  *  ls 
and  Vp||  -  1.004C,  the  lowbltron  mode  becomes 
dominated  over  the  electron  cyclotron  mode  and  the 
electron  beam  is  able  to  convert  25Z  of  Its 
kinetic  energy  into  radiation. 


Efficiency  Enhancement 


that  the  efficiency  can  be  enhanced  from  24Z  to 
42Z.  Por  competition,  the  time  history  of  elec¬ 
tron  phase  space  Is  shown  In  Pig.  4.  at  the  time 
#1  saturation  (*ct  ■  145)  for  the  case  of  varying 
aolenoldal  field,  the  electrons  are  mostly  in 
their  decelerating  phase  (Pig.  4c),  whereas  the 
electrons  are  mostly  In  their  accelerating  phase 
(Pig.  4b)  without  varying  B0. 

Summary 

In  conclusion,  by  choosing  the  parameters  such 
that  Vpj,  c  and  Bw  ■  B0  a  lowbltron  la  capable  of 
attaining  high  efficiency  of  the  order  of  2SZ. 
Furthermore,  the  efficiency  can  be  Increased  to 
42Z  by  properly  profiling  the  solenoldal  magnetic 
field. 
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A  new  scheme  of  enhancing  the  efficiency  and  gain  of  a  Doppler-shifted  electron-cyclotron  autoreso¬ 
nance  maser,  with  use  of  a  transverse  magnetic  wiggler  to  inhibit  stimulated  absorption,  is  proposed. 
Computer  simulations  with  use  of  a  if  -dimensional,  electromagnetic,  relativistic  particle  code  show 
significant  improvement  in  the  efficiency,  and  gain  may  be  attained  with  this  scheme. 

PACS  numbers:  42.J2.+* 


Advances  in  high-power,  relativistic  electron  beam 
technology  over  the  past  decade  have  made  the  electron- 
cyclotron  maser1  -2  an  increasingly  attractive  source  of 
tunable,  coherent  radiation  in  the  millimeter  and  submil¬ 
limeter  wavelength  ranges.  The  radiation  is  generated 
by  relativistic  electrons,  which,  as  a  result  of  the  depen¬ 
dence  of  electron-cyclotron  frequency  on  the  relativistic 
mass,  tend  to  bunch  together  as  they  gyrate  around  a 
solenoidal  magnetic  field,  B»3  The  resonant  condition 
for  electron-wave  interaction  is  given  approximately  by 

(or  “*ri!+  fla/y,  (1) 

where  "j  is  the  drift  velocity  of  the  beam  along  Bo, 
y“(l  —  c  2/c 2 ) l/2,  fto/y  is  the  relativistic  electron-cyclo¬ 
tron  frequency,  and  is  the  phase  velocity  of  the 

resonant  wave.  An  output  frequency  much  higher  than 
the  electron-cyclotron  frequency  may  be  obtained  by  use 
of  highly  relativistic  beam  drift  velocity.'1 

A  property  of  electron  interaction  with  the  wave  is 
that  its  velocity,  r(,  is  a  function  of  its  energy  and  initial 
conditions  only,  regardless  of  the  electromagnetic  (EM) 
wave  amplitude2: 

0i -ll  “(1  “0ph/3,o)y</r)/0ph:  (2) 

where  fii^ci/c,  fiph  "‘ta/kc,  yo*“y(to)>  and  fino^fiito). 
This  follows  from  the  equations  for  the  electron  energy 
and  momentum,  dfydt  “  —  ev  E  and  dpjdt  - -e(v 
•E)/fp *.  Equation  (2)  shows  that  an  electron  can  lose 
all  its  kinetic  energy,  both  transverse  and  longitudinal,  to 
the  wave  if  its  initial  drift  velocity  satisfies  fiw 
mm(ro~  I  )/ro/V  In  general,  electrons  confined  to  ellip¬ 
soidal  trajectories  in  fix"  fit  phase  space 

-((I  “0ph0i)/K>(l  -0ph0iio))J.  (3) 

suffer  maximum  energy  loss  when  they  have  given  up  all 
their  transverse  momentum.  For  efficient  energy  conver¬ 
sion,  an  electron  has  to  be  able  to  maintain  the  resonance 
condition  for  a  sufficiently  long  time.  The  phase  shift  of 
an  electron  relative  to  the  wave  may  be  expressed  as 
o(r)  “(6*  +  Sj).  where  8k  ™(o>  —  <o,)t  is  the  kinetic  shift 
and  8j  "  —  Uac,+  rigAO/fllf  is  the  dynamic  shift  due 


to  electron-wave  interaction.  With  use  of  Eqs.  (I)  and 
(2),  Sj  may  be  expressed  as 

64  ■  [a>,(Ay/yo)Hl  ~  (co/ajrPph)]t.  (4) 

The  dynamical  shift  vanishes  if  fifa-w/to,,  the  change  in 
electron-cyclotron  frequency  being  exactly  cancelled  by 
the  change  in  Dopper  shift  due  to  electron  recoil.  There¬ 
fore,  if  to^to,  is  initially  satisfied,  an  electron  will 
remain  in  phase  with  the  EM  wave  for  all  times  if 
fip h"l.  This  autoresonant  regime  holds  out  the  prospect 
of  high  efficiency. 5  6 

However,  depending  on  the  initial  phase  of  electrons 
relative  to  the  wave,  the  autoresonant  condition  can  also 
lead  to  acceleration  of  electrons.  Equations  (2)  and  (4) 
show  that  an  electron  can  be  accelerated  indefinitely. 
The  pitch  factor  of  the  beam  is  given  by  0^/0, 
"{2/(y- 1 )] 1/2  (for  0Ph™l).  so  that  at  large  y,  the 
beam  momentum  is  predominantly  along  Bo-  This  ac¬ 
celeration  mechanism  has  attracted  atten’ion  to  its  po¬ 
tential  application  as  a  high-energy  particle  accelerator.7 

The  electron-energy  loss  and  gain  are  the  classical 
correspondents  of  stimulated  emission  and  absorption. 
To  encourage  positive  gain  in  the  cyclotron  autoreso¬ 
nance  maser  (CARM),  a  small  initial  mismatch  in  fre¬ 
quency,  to  >  a >„  with  0ph  >  <o/ai ,  is  desirable  in  order  to 
favor  stimulated  emission  over  absorption.2  An  exam¬ 
ination  of  the  equations  for  phase  shifts  above  shows  that 
electrons  that  give  up  energy  to  the  wave  will  suffer 
smaller  phase  shifts  than  those  that  are  accelerated  by  it. 
with  these  conditions.  Nevertheless,  as  the  wave  grows 
in  a  conventional  CARM,  a  sizable  fraction  of  the  elec¬ 
trons  will  be  accelerated  by  it  and,  since  wave  energy  is 
expended  in  accelerating  these  electrons,  they  depress 
both  the  efficiency  and  gain  of  the  CARM. 

In  this  paper  we  propose  a  method  of  curtailing  the 
autoresonant  acceleration  of  electrons.  We  show  that  bv 
imposing  a  transverse,  helical  “wiggler”  magnetic  held 
onto  a  conventional  CARM,  electrons  accelerated  by  the 
wave  can  be  made  to  go  off  resonance  with  the  wave, 
and,  hence,  terminate  the  accleration  process,  while  elec¬ 
trons  that  lose  energy  to  the  wave  will  remain  in  reso¬ 
nance  with  it.  Computer  simulations  demonstrate  that 
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both  efficiency  and  growth  rate  can  be  significantly 
enhanced  as  a  result. 

The  detuning  mechanism  arises  from  the  behavior  of 
electrons  moving  in  the  combined  wiggler  and  uniform 
axial  magnetic  field  modeled  by 

+  licos(/t»r  )+ysin(A»r)l, 

where  B »  is  the  Higgler  field  strength  and  A»  its  pilch. 
This  wiggler  field  is  an  idealization  of  that  produced  by  a 
bifilar  curreni  winding  and  currently  in  use  in  Tree- 
electron  laser  experiments.'  Equilibrium  electron 
motion  in  this  magnetostatic  field  is  approximately  de¬ 
scribed  by 

cos(A».z)  +  iJ.cos(n<-r  +  p), 

«>  *i «  sin(A.r)4i  Asin(n,r +p),  (5) 

i'.-  “i'b  ~pcos(A,2  —  nri  +  p), 

where  r.  -r,(Bm/B0)/(l -kwrt/n, ),  p -ivrA/ft.  nf 
■•no /y,  0  is  an  initial  phase,  and  'p/in«c I  has  been  as¬ 
sumed.  The  electron  beam  is  injected  with  a  finite  trans¬ 
verse  velocity,  r  A,  into  the  wiggler  region.  The  effect  of 
B.  is  to  induce  an  additional  mode  of  gyration  about  2 
with  a  frequency  n  =» and  an  oscillation  in  the  axial 
velocity  r,  about  its  mean  value  n,  at  the  beat  frequency 

From  Eqs.  (2),  (3),  and  (5),  we  can  deduce  approxi¬ 
mately  the  variation  of  p  with  i  ».  This  is  shown  qualita¬ 
tively  in  Fig  I.  A  singularity  occurs  at  gyrorcsonance, 
Anr(mne,  in  the  neighborhood  of  which  Eqs.  (5)  no 
longer  give  an  accurate  description  of  electron  motion. 
On  either  side  of  the  singularity,  p  varies  monotonicaily 
with  iv.  branch  Af  corresponds  to  the  case  flf  >  A„r»  and 
branch  f  to  flf  <  A,  tv  Figure  I  is  also  qualitatively 
representative  of  the  variation  of  p  with  y,  with  the 
transformations  y-*  I  as  fa  —  0  and  y— *  «  as  —  I. 
Assumed  in  this  qualitative  description  of  p  is  that  the 
wiggler  is  sufficiently  weak  so  that  the  mean  B  is  still 
essentially  given  by  Eq.  (2).  This  condition  is  necessary 
if  we  are  to  use  the  cyclotron  autoresonance  effect.  Also 


FIG  I.  Variation  of  amplitude  of  axial  oscillation  p/r  with 
mean  velocity.  0, 


assumed  is  that  the  electron  motion  can  be  approximate¬ 
ly  described  by  Eq.  (5)  even  in  the  presence  of  the  EM 
field  (see  a  discussion  of  Fig.  3  below). 

The  variation  of  p  with  r«  (and  y)  depicted  in  Fig  I 
can  be  exploited  to  provide  a  discriminate  detuning 
mechanism  in  the  CARM  Consider  a  CARM  with  a 
helical  magnetic  field  such  that  initially  nf>A„ii  and 
the  electron  beam  is  described  by,  for  example,  point  A 
on  curve  X.  Then  electrons  accelerated  by  the  wave  will 
move  to  the  right  of  A,  gaining  a  larger  oscillating  com¬ 
ponent  in  its  axial  velocity.  If  the  slope  at  A  is 
sufficiently  steep,  the  increasing  axial  oscillations  will 
quickly  desynchronize  these  electrons  with  the  wave.  On 
the  other  hand,  electrons  losing  energy  to  the  wave  will 
move  to  the  left  of  A,  reducing  their  axial  oscillations 
and.  hence,  becoming  more  in  resonance  with  the  wave. 
In  short,  accelerated  electrons  become  more  “noisy"  and 
quickly  move  out  of  resonance  with  the  wave,  while  the 
opposite  process  occurs  for  those  electrons  that  are  driv¬ 
ing  the  wave. 

To  illustrate  the  mechanism  behind  this  enhancement 
scheme  and  its  potential,  we  present  below  results  from 
computer  simulations  of  a  CARM  with  helical  magnetic 
field  using  a  1  } -dimensional,  EM,  relativistic  particle 
code  with  periodic  boundary  conditions.9  The  electrons 
are  initially  monoenergetic  with  uniform  axial  velocity, 
i'i.  and  vx,  and  are  uniformly  distributed  over  the  trans¬ 
verse  phase  angle,  e“,lan-l(i-,/i'>),  —  * <6<n.  The 
helical  field,  is  switched  on  adiabalically  to  let  elec¬ 
trons  find  their  equilibrium  orbits.  The  EM  wave  is  then 
allowed  to  grow  from  noise,  with  B»  kept  constant.  Only 
one  A  mode  cf  the  EM  wave  is  kept  in  a  given  simula¬ 
tion.  Results  are  presented  below  for  the  following  pa¬ 
rameters  and  initial  conditions: 

yo“3,  B»/Bq“0.09) ,  no/yo^»t'iio"*2.69, 

fto/yotup”2.33,  0xo"r<T\ 

where  a>r  ~(4jrnee2/m)]/2.  The  values  correspond  to  a 
CARM  operating  along  curve  X  in  Fig.  1.  Comparisons 
are  made  with  simulations  with  identical  initial  condi¬ 
tions  but  with  B»  “0.  As  BJB o<K  1,  and  the  initial  con¬ 
ditions  are  removed  from  gyroresonance,  B»  initially 
presents  only  a  small  perturbation  to  conventional 
CARM  electron  motion. 

Figure  2  compares  efficiency  and  growth  rate  with  and 
without  B„  as  a  function  of  fa,  the  phase  velocity  of  the 
excited  EM  wave.  The  efficiency  is  defined  as  r?“(y<, 
-y)/(ro_  1)-  where  y  is  the  value  of  y  averaged  over 
ail  electrons  at  saturation.  Substantial  improvement  in 
gain  and  growth  rate  is  evident  over  a  range  of  fa 
slightly  greater  than  1.  The  peak  efficiency  is  increased 
by  about  a  factor  of  3  with  the  inclusion  of  B ».  from  1 1^ 
to  3291.  The  effectiveness  of  B.  in  suppressing  au- 
toresonant  acceleration  can  be  seen  from  Fig  3.  which 
shows  phase  plots  of  kinetic  energy  {i  “(y~  1 )  vs  p  for 
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FIG.  2.  Comparison  of  (a)  efficiency  and  (b)  growth  rate, 
with  (filled  circles)  and  without  (open  circles)  8.,  as  a  func¬ 
tion  of  the  phase  velocity  fa  —  ai/kc. 
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FIG.  3.  Phase  plots  of  (y —  ' )  vs  0,  at  saturation  for 
Pph~  I  002  (a)  with  8.  and  (b)  without  8.. 


the  case  PPh  =  1  002  with  and  without  Bw  (the  second  set 
of  data  points  from  the  left-hand  side  in  Fig.  2(a))  at 
saturation.  Without  B«.  electrons  are  distributed  along 
the  h  i  bola  defined  by  Eq.  (2).  Initially,  all  electrons 
lie  i  point  (£*  ■■2,  0"0,o)  on  this  curve.  As  the  EM 

wave  .rows,  electrons  driving  the  wave  move  down  the 
curve,  while  :hose  accelerated  by  the  wave  move  up  the 
curve.  The  lower  boundary  of  this  curve  is  defined  by 
the  value  of  0,  at  which  an  electron  loses  all  transverse 
momentum.  For  this  example,  this  occurs  at  0q=*O.78, 
&a«0.6,  from  Eq.  (3).  At  the  upper  end  of  this  curve, 
electrons  are  accelerated  up  to  **  4.3.  Acceleration 
eventually  stops  when  phase  slippage  between  electron 
and  wave  exceeds  ±  x.  For  the  case  with  B„,  electrons 
are  initially  described  by  a  line  joining  the  points  ({*  “2, 
01-0.82)  and  ($*—2,  0,  "0.92).  This  is  because  the 
monoenergetic  electrons  are  oscillating  about  their  mean 
axial  velocity  such  that  A0»— 2ji.  As  electrons  lose  ener¬ 
gy  to  the  wave,  n  becomes  smaller,  as  discussed  above, 
and  eventually  vanishes  when  r  x  *0,  hence  the  V  shape 
of  the  phase  space  enclosing  the  electrons.  The  width  of 
this  V-shaped  region  as  a  function  of  r  agrees  closely 
with  theoretical  prediction  based  on  Eqs.  (2),  (3),  and 
(3),  showing  that  these  equations  do  describe  the  elec¬ 
tron  motion  adequately  in  this  example.  Maximum  en¬ 
ergy  loss  for  an  electron  again  corresponds  approximate¬ 
ly  to  *0.6  and  0,4*0.78.  However,  electron  accelera¬ 


tion  has  been  drastically  curtailed.  Maximum  electron 
energies  are  only  about  '"3.0.  less  than  half  that 
gained  without  Bw. 

It  is  possible  to  obtain  efficiencies  higher  than  the  1 1  rc 
here  for  conventional  CARM's  (uniform  solcnoidal  field 
with  no  wiggler)  by  use  of  other  operating  parameters. 
Efficiencies  in  excess  of  20%  have  been  shown  to  be  pos¬ 
sible  in  semianalytic  calculations  with  initial  conditions 
that  yield  a  single-panicle  efficiency  of  100%. 4  The 
single-particle  efficiency  (the  maximum  efficiency  an 
electron  can  achieve)  may  be  obtained  from  Eq.  (3)  by 
putting  0X  “0  (the  figure  for  the  above  example  is  70%). 
And  in  3D  wave-guide  simulations,  efficiencies  of  about 
30%  have  been  obtained. 5  Our  intention  in  this  Letter  is 
only  to  point  out  that  the  use  of  a  wiggler  can 
significantly  enhance  the  CARM.  The  maximum 
enhancement  possible  for  an  optimally  designed  conven¬ 
tional  CARM  will  be  explored  in  later  studies.  Aiso  to 
be  determined  are  the  scalings  of  Bw  and  kw  with  the 
operating  parameters  of  a  CARM  (Bo.  Fo.  0«o)  for  op¬ 
timum  enhancement.  Their  values  here  were  chosen  by 
a  process  of  trial  and  error.  A  judicious  choice  of  8* 
and  kw  is  necessary  because  if.  for  example.  £„  is  too 
large,  the  resulting  large  axial  oscillations  in  the  elec¬ 
trons’  motion  will  prevent  the  resonance  condition  (Eq. 
(1))  from  being  satisfied  by  any  electron,  and,  hence,  no 
gain  will  result,  while  if  Bw  is  loo  small,  the  enhancc- 
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mcnt  will  be  unattractive.  As  a  rough  measure  of  (he 
sensitivity  to  B%  in  this  example,  decreasing  5.  by  50r< 
reduces  the  efficiency  of  (he  A  mode  corresponding  to  the 
peak  efficiency  of  325  to  195,  while  increasing  B,  by 
255  and  505  reduces  the  efficiency  to  255  and  less  than 
15.  respectively. 

We  note  that  this  scheme  is  different  from  previous 
proposals  for  enhancement  which  are  based  on  tapering 
the  solenoidal  field.  One  merit  of  it  is  that  it  does  not  re¬ 
quire  the  development  of  any  new  technology.  The  heli¬ 
cal  magnetic  field  modeled  is  already  in  use  as  the 
wiggler  in  many  free-electron-laser  (FEL)  experiments. 
But  unlike  FEL's,  CAR  M’s  require  the  electron  beam  to 
have  an  initial  transverse  momentum,  which  may  be 
effected  by  passing  the  beam  through  a  wiggler  or  along 
a  nonuniform  solenoidal  field.  One  problem  has  been 
that  the  high  uniformity  in  transverse  momentum  and 
low  noise  level  desirable  for  optimum  operation  have 
proved  difficult  to  achieve.  However,  the  enhancement 
in  gain  and  efficiency  that  may  be  possible  by  use  of  a 
wiggler  may  help  to  alleviate  the  problem  of  achieving 
high  beam  quality  and  increase  the  attractiveness  of  the 
CARM  over  the  FEL  which  typically  has  an  efficiency  of 
only  a  few  percent  (with  a  uniform  wiggler). 

This  scheme  may  also  have  potential  application  in  a 
cyclotron  autoresonance  particle  accelerator. 7  It  may  be 
seen  from  Fig.  1  that  selective  suppression  of  au- 
loresonant  deceleration  can  be  achieved  by  the  same 


means  as  autoresonant  acceleration.  By  choosing  the  ini¬ 
tial  conditions  so  that  A„r*o>  fWyo  (curve  )'),  electrons 
losing  energy  to  the  wave  can  be  made  to  go  off  reso¬ 
nance  with  the  wave  rapidly,  hence  terminating  the  de¬ 
celeration. 
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ABSTRACT 


Theoretical  analysis  and  computer  simulations  have  been  carried  out 
to  examine  the  stability  and  tunability  of  a  CARM  (cyclotron  auto  resonance 
maser)  amplifier.  It  is  shown  that  frequency  tunability  over  one  octave 
can  be  achieved  by  magnetic  field  tuning  within  the  stabilllty  range. 
Tunability  can  be  further  enhanced  if  the  waveguide  wall  is  made  lossy 
to  increase  the  threshold  of  the  absolute  instability.  However,  beam 
momentum  spread  is  found  to  have  a  strong  deteriorating  effect  on  the 
tunability  and  gain. 


KWS 


I.  Introduction 


During  the  last  two  decades,  there  is  intense  interest  in  developing 
fast-wave  devices  based  on  the  mechanism  of  electron  cyclotron  resonance 
maser.  The  resonant  condition  between  an  electron  and  an  electromagnetic  wave 
propagating  in  uniform  magnetic  field  is  o>  ■  ftc  +  kzvz  where  m  is  the  wave 
frequency  kzvz  is  the  Doppler  shift,  and  is  the  relativistic  electron 
cyclotron  frequency.  Gyrotron  [1]  operates  in  the  regime  where  »  kzvz. 
The  efficiency  of  a  relativistic  gyrotron  is  lower  than  a  weakly  relativistic 
gyrotron  because  in  the  former  case  a  small  percentage  change  in  the 
electron  energy  will  result  in  a  large  detuning  of  the  resonant  condition. 

In  the  case  of  a  gyro-TWT  (gyrotron  travelling  wave  tube  amplifier),  the 
operating  frequency  is  only  slightly  above  the  cutoff  frequency  of  the 
waveguide  and  hence  is  very  close  to  the  oscillation  frequency  of  the  absolute 
Instability.  This  makes  it  difficult  to  find  a  scheme  to  suppress  the 
absolute  instability  without  affecting  the  amplified  wave. 

Most  of  these  problems  can  be  alleviated  by  going  into  the  regime  where 
kzvz  dominates  over  [3,4,].  Substituting  Vph  (  =  <*>/kz)  for  kz,  the 
resonance  condition  becomes  w  =  nc/(l-  vz/vph).  Thus  w  can  be  substantially 
higher  than  flc  if  vz  is  close  to  vph»  It  has  also  been  shown 
that  under  the  condition  vph  <  c 2/vz,  the  increase  in  fie  can  be  compensated 
by  the  decrease  in  k2vz  such  that  the  resonance  condition  is  maintained  over  a 
long  interation  time.  This  is  the  auto  resonance  conditon  exploited  by  the 
CARM.  In  a  CARM  amplifier,  the  operating  frequency  is  well  above  the  cutoff 
frequency.  When  oscillation  does  take  place  near  the  cutoff  frequency  due  to 
an  absolute  Instability,  it  is  relatively  easy  to  selectively  supress  the 
absoulute  instability.  In  addition,  the  oscillation  can  be  suppressed  either 
by  injecting  a  signal  with  intensity  well  above  the  noise  level  [5]  or  by 


resiatively  coating  the  waveguide  wall. 

In  addition  to  the  advantages  of  high  efficiency  and  high  frequency 
(relative  to  Ac),  the  frequency  of  a  CARM  amplifier  is  tunable  by  means  of 
varying  either  the  electron  beam  energy  or  magnetic  field.  In  this  paper  the 
latter  scheme  is  investigated  through  theoretical  analysis  and  computer 
simulations.  Instead  of  covering  a  broad  data  base,  we  will  illustrate  the 
general  properties  through  a  particular  example  given  in  Table  I.  The  results 
show  a  frequency  tunability  of  more  than  100%  over  a  magnetic  field  range  of 
40%  (within  this  range  the  amplifier  remains  absolutely  stable).  The 
effects  of  beam  momentum  spread  and  waveguide  wall  resistivity  are  also 
addressed.  The  linear  theory  is  presented  in  Section  II,  which  is  found  to  be 
in  good  agreement  with  computer  simulations  given  in  Section  III. 


II.  Linear  Theor 


Brigg's  criteria  for  absolute  and  convective  Instabilities  [6]  and 
previous  linear  theories  of  the  gyro-TWT  [2,7,8]  have  been  applied  to  the  case 
under  study. 

Fig.  la  shows  the  intersection  of  the  waveguide  mode  and  beam-wave 
resonance  line  at  a  grazing  angle.  This  is  typical  of  the  gyro-TWT  operation. 
For  the  CARM,  it  is  preferable  to  operate  at  a  higher  magnetic  field  so  that 
the  upper  Intersecting  point  has  a  phase  velocity  near  the  speed  of  light. 

This  will  result  in  auto  resonance  and  therefore  the  extraction  of  beam 
streaming  energy.  When  the  magnetic  field  is  raised,  however,  the  lower 
Intersecting  point  will  move  to  the  left  and  eventually  an  absolute 
instability  (oscillation)  sets  in  as  a  result  of  backward  wave  excitation. 
Thus,  for  the  stable  operation  of  the  CARM,  there  exists  an  upper  limit  on 
the  magnetic  field,  which  we  label  as  Bosc. 

When  the  magnetic  field  is  above  the  grazing  magnetic  field  Bg  (for  which 
the  grazing  Intersection  in  Fig.  la  is  realized),  the  approximate  analytical 
criteria  for  the  prediction  of  absolute  instability  derived  in  Ref.  2  is  no 
longer  valid.  Here,  we  resort  to  the  original  criteria  in  Ref.  6  for  the 
determination  of  the  threshold  of  the  absolute  instability  by  numerically 
tracing  the  behavior  of  the  roots  of  the  dispersion  relation  [7]  in  the  o>— k.z 
space.  Fig.  lb  plots  Bosc/Bg  as  a  function  of  a  wall  loss  parameter  5,  defined 
as  the  ratio  of  the  wall  skin  depth  to  the  wall  radius.  It  is  seen  that  there 
is  a  40-602  range  of  stable  magnetic  field,  depending  on  the  wall  loss.  We 
note  that  6  *  0.01  represents  an  extremely  lossy  wall,  which  can  be  realized 
only  for  waveguide  of  very  small  radius.  A  magnetic  field  tuning  range  of 
402,  however,  can  lead  to  a  frequency  tuning  range  over  100%,  as  is  shown  in 


Also  shown  in  Pig.  2  is  the  linear  gain  in  one  cutoff  wavelength  (2irc/wc) 
of  the  waveguide  for  four  values  of  axial  momentum  spread  Apz/p2Q,  where  Apz 
is  the  standard  deviation  of  pz  from  the  average  value  pzQ,  assuming  a 
Gaussian  type  distribution.  The  gain  is  more  sensitive  to  the  momentum  spread 
at  higher  magnetic  field  (hence  higher  frequency).  This  is  because  velocity 
spread  results  in  broader  spread  in  Doppler  shifts  at  larger  propagation 
constant  kz.  Because  of  its  sensitivity  to  the  velocity  spread,  the  CARM 
amplifier  requires  a  very  good  quality  beam  in  order  to  operate  in  a  highly 
Doppler  shifted  regime. 

If  wall  loss  is  enhanced  to  achieve  higher  tunability  (Fig.  lb),  a 
fraction  of  the  extracted  beam  energy  will  be  dissipated  on  the  wall.  Fig.  3 
shows  that,  in  the  preferred  regime  of  CARM  operation  (high  end  of  the  stable 
magnetic  field  range),  the  Ohmic  power  loss  is  small  as  compared  with  the 


output  wave  power 
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III.  Computer  Simulation  Results 

Particle  simulation  is  a  very  powerful  tool  to  model  various  free  electron 
devices  [9].  In  our  simulation  code  the  transverse  profile  of  the  electro¬ 
magnetic  fields  is  fixed  (TEqi  mode)  whereas  the  longitudinal  profile  and 
temporal  evolution  are  self-consistently  determined  by  Integrating  the 
Maxwell's  equations  using  the  perturbed  beam  current  as  the  source. 
Longitudinal  dependence  of  the  fields  is  decomposed  into  Fourier  components 
to  provide  resolution  of  the  simultaneous  presence  of  different  modes.  This 
is  an  essential  feature  for  analyzing  the  absolute  instability  in  a  CARM 
amplifier.  Electrons  are  injected  at  the  entrance  end  with  prescribed 
parameters  and  diagnosed  at  the  output  end.  Their  trajectories  in  the 
Interaction  region  are  followed  through  the  relativistic  equation  of  motion 
and  the  Lorentz's  forces  consist  of  the  self-consistent  as  well  as  externally 
applied  fields.  Typically,  thousands  of  electrons  are  used  so  that  beam 
momentum  spread  is  well  represented.  The  fact  that  electrons  experience 
inhomogeneous  electromagnetic  fields  is  included  In  the  code  by  allowing 
electron  motion  in  three  spatial  dimensions.  The  wall  resistivity  is  modelled 
by  imposing  a  damping  term  in  the  Maxwell's  equation  with  a  damping  rate  g/o^ 

■  5/rw  where  6  and  rw  are  respectively  the  skin  depth  and  waveguide  radius. 
Using  the  parameters  in  Table  I,  a  series  of  simulations  have  been  carried  out 
to  study  the  stability  and  tunability  of  a  CARM  amplifier. 

Fig.  4a  shows  the  steady  state  (<oct  -  200)  spatial  profile  of  the  wave 
field  for  a  run  with  Apz/pzo  ■  0  and  B  -  1.15Bg.  The  drive  power  Is  3.75  MW 
at  a  frequency  of  2.663  uc.  Immediately  after  lanunchlng,  the  wave  shows  no 
appreciable  growth  because  of  the  launching  loss.  There  is  then  a  region  of 
exponential  growth  giving  a  linear  growth  rate  of  approximately  3dB  per  cutoff 
wavelength,  in  very  good  agreement  with  the  analytical  prediction.  Saturation 


is  reached  at  the  output  end  where  the  efficiency  is  approximately  20%.  Fig.  4b 


shows  the  electron  kinetic  energy  versus  its  phase  ($  ■  o>Qt  -  kzz-tan“l— ) 

Py 

relative  to  the  signal  wave.  The  saturation  is  caused  by  phase  trapping  of 
the  electrons  as  can  be  clearly  seen  in  the  phase  space  plot  that  at  the  exit 
plane  where  an  electron  bunch  is  formed  in  the  decelerating  phase  and  it  Is 
rotating  in  phase  space.  Two  other  simulations  with  B-Bg  and  1.3  Bg  (which 
are  still  in  the  stable  regime  according  to  the  theoretical  prediction)  have 
been  performed.  The  efficiencies  are,  respectively,  16%  and  23%.  The  linear 
gains  for  all  three  cases  are  plotted  in  Fig.  3  as  solid  dots.  To  investigate 
the  effects  of  beam  momentum  spread  on  the  linear  gain,  the  three  cases  are 
rerun  with  Apz/pzQ  increased  to  0.01.  The  results  are  also  shown  in  Fig.  3. 
The  spread  substantially  reduces  the  linear  gain,  especially  in  the  high 
magnetic  field  regime  of  interest  to  GARM.  The  devices  become  unattractrlve 
due  to  low  gain  and  low  efficiency  if  beam  momentum  spread  is  greater  than  3%. 

I 

In  the  case  of  B»Bg  (gyro-TWT  regime),  the  phase  velocity  of  the 
signal  wave  is  too  far  above  the  speed  of  light  to  satisfy  the  auto  resonance 
condition.  Hence,  only  the  transverse  beam  energy  is  converted  into 
radiation,  as  can  be  seen  in  Fig.  5a  where  the  time  evolution  of  beam  trans¬ 
verse  and  longitudinal  energies  at  the  output  end  is  shown.  As  the  magnetic 
field  is  raised  to  1.3  Bg  the  phase  velocity  of  the  signal  wave  becomes  close 
enough  to  the  speed  of  light  to  satisfy  the  auto  resonance  condition.  This 
results  in  the  conversion  of  both  beam  transverse  and  longitudinal  energies 
into  radiation  which  is  clearly  exhibited  in  Fig.  5b.  The  saturated 
efficiency  therefore  increases  from  16%  in  the  gyro-TWT  regime  to  23%  in  the 
CARM  regime.  According  to  the  theory  in  Sec.  II  (Fig.  16),  further  increasing 
the  magnetic  field  to  B>  1.4Bg  will  bring  in  the  absolute  instability.  A  run 

with  B  “  1.5  Bg  has  been  carried  out  to  confirm  this  prediction.  The  wave 
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spatial  profiles  at  two  instants  are  shown  in  Fig.  6a  (ti>ct«100)  and  Fig.  6b 
(«ct-200).  As  the  absolute  instability  grows  into  long  wavelength 
oscillation,  it  modulates  the  amplified  wave  and  results  in  a  reduction  in 
output  power  at  the  desired  frequency.  One  method  of  suppressing  oscillations 
is  by  resistive  coating  of  the  waveguide  wall.  This  scheme  allows  the  device 
to  operate  at  a  higher  magnetic  field  and  thereby  enhances  the  efficiency  and 
Increases  the  tunable  range.  Zn  order  to  investigate  these  effects,  a  run 
with  B  -  1.5  Bg  and  6  -  0.08  rw  is  performed  and  the  spatial  distribution  of 
the  wave  field  at  o»ct  -  200  is  displayed  in  Fig.  6c.  The  wall  loss 
substantially  but  not  completely  suppressed  the  long  wavelength  oscillation 
with  negligible  effect  on  the  amplified  wave.  The  output  power  in  Fig.  be  is 
mainly  at  the  input  frequency.  However,  the  required  resistivity  for  absolute 
instability  stabilization  is  almost  one  order  of  magnitude  higher  in 
simulation  runs  than  predicted  by  theory. 

In  summary,  theoretical  analysis  and  computer  simulations  have  been 
carried  out  to  investigate  the  stability  and  tunability  of  a  CARM  amplifier. 
The  results  show  that  more  than  100%  tunability  can  be  achieved  by  magnetic 
field  tuning  within  the  stability  region.  Beam  quality  is  of  critical 
importance  to  the  performance  of  the  CARM  amplifier. 
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Figure  Captions: 

Fig.  1  (a)  Grazing  Intersection  of  the  waveguide  mode  and  the  beam-wave 

resonance  line,  (b)  Bosc/Bg  versus  the  wall  loss  parameter  5. 
Parameters  In  Table  I  were  used  for  Fig.  lb  and  subsequent  figures. 

Fig.  2  Center  frequency  (dashed  line)  and  linear  gain  per  cutoff  wavelength 
(solid  lines)  versus  the  operating  magnetic  field.  5  -  0  is  used. 
Gains  for  5  -  0.005  and  0.01  are  only  slightly  different  (less  than 

5Z  lower).  Simulated  linear  gain  (dots)  is  in  good  agreement  with 
theory. 

Fig.  3  Power  lost  on  the  wall  as  a  function  of  the  operating  magnetic  field 
The  corrersponding  operating  frequency  is  shown  in  Fig.  2. 

Fig.  4  Simulation  results  for  B  «  1.15  Bg  (a)  spatial  profile  of  the 

wave  field  at  «ct  ~  200.  (b)  electron  kinetic  energy  versus  its 

phase  relative  to  the  input  wave. 

Fig.  5  Simulated  evolution  of  electron  beam  axial  and  transverse  energies 
at  the  output  enld  for  a)  B  ■  Bg,  (b)  B  »  1.3  Bg. 

Fig.  6  Simulated  spatial  profiles  of  the  wave  field  for  B  ■  1.5  Bg. 

(a)  wct  =.  100.  (5  =*  0)  (b)  met  -  200.  (6-0)  (c)  <uct  -  200. 
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Introduction; 

This  report  describes  the  work  done  by  Varian  Associates  on  the  study 
of  the  Cyclotron  Auto  Resonance  Maser  (CARM)  for  the  contract  with  U.C.L.A. 
during  the  fiscal  year  1986.  The  approach  taken  was  to  study  the  CARM 
oscillator  from  the  point  of  view  of  designing  a  specific  experiment, 
schematically  illustrated  in  Figure  1.  The  design  parameters  chosen  would 
allow  the  use  of  an  existing  500  kV,  200  amp  Magnetron  Injection  Gun 
design  (with  some  modification)  and  a  28  GHz  gyrotron  magnet.  This  could 
allow  a  cost  effective  experiment  to  be  performed  which  would  be  sufficient 
to  demonstrate  the  CARM  as  a  viable  device  and  verify  the  basic  physics. 

The  rest  of  this  report  is  organized  as  follows:  Section  I  describes 
how  the  basic  design  parameters  where  chosen.  In  Section  II,  prediction  of 
saturated  output  power  and  efficiencies  are  presented  indicating  regions  of 
hard  and  soft  excitation.  Section  III  gives  preliminary  designs  for  the 
magnetron  Injection  Gun  and  Bragg  Reflector.  Finally,  Section  IV  discusses 
the  future  work  required  to  improve  and  fully  verify  the  final  design. 


Section  I.  Design  Considerations 


The  beam  voltage  and  current  were  chosen  to  be  500  kV  and  200  amps 
respectively  in  order  to  make  use  of  an  existing  gyroklystron  gun  design. 
The  cutoff  frequency  of  the  resonator  section  was  kept  at  10  GHz  in  order 
to  keep  the  final  beam  guiding  center  the  same  as  in  the  gyroklystron.  This 

also  required  staying  with  the  TEq^  mode.  The  cyclotron  frequency  and 
operating  frequency  were  determined  from  the  conditions  for  cyclotron 
resonance  and  ideal  cavity  resonance  according  to  the  following  formula: 


co 


1 


^  vz/vph 


i 


a 


“c  f- 1  -  ^/’ph^ 

where  Vp^  and  vz  are  the  phase  velocity  and  axial  beam  velocity.  Previous 
work  done  at  U.C.L.A.  served  as  a  guide  in  choosing  a  phase  velocity  of 
2.69%  above  the  speed  of  light  in  order  to  allow  a  reasonable  compromise 
between  gain  and  efficiency.  This  implied  an  operating  frequency  of  44  GHz. 
The  transverse  velocity\c  was  chosen  slightly  greater  than  1 lY  (enhanced 
gain)  resulting  in  a  ratio  of  perpendiclar  to  parallel  velocity  (e< )  of  .8. 
This  resulted  in  a  ratio  of  cyclotron  frequency  to  operating  frequency  of 
M/3  translating  into  a  required  axial  magnetic  field  of  10300  gauss.  The 
cavity  interaction  length  was  determined  by  calculating  oscillation 
threshold  currents.  The  calculation  was  done  assuming  the  fields  were  of 
an  ideal  cavity  resonator  whose  length  was  N  half  wavelength  long.  This 
was  considered  for  now  to  be  a  reasonable  approximation  to  the  true  open 
resonator  which  would  be  formed  from  Bragg  reflectors.  The  cavity  Q  was 


assumed  to  be  500.  The  results  shown  in  Figure  2  indicated  that  lengths 
greater  than  50  half  wavelengths  would  be  required  in  order  to  allow  start 
oscillation  currents  of  less  than  200  amps.  An  optimal  design  with  a 
length  of  80  half  wavelengths  and  a  Q  of  650  was  finally  chosen  although 
much  of  the  following  analysis  was  done  with  a  cavity  length  of  70  half 
wavelengths  an<i  a  Q  of  600.  A  summary  of  the  basic  design  parameters  is 
shown  in  Table  I. 

Section  II.  Predictions  of  Saturated  Output  Power  and  Efficiency 

The  efficiency  analysis  was  performed  using  a  simplified  semi- 
analytical  model  which  assumed  that  the  saturated  steady  state  longitudinal 
mode  structure  was  determined  by  geometry  only  and  not  greatly  perturbed  by 
the  presence  of  the  beam.  A  single  equation  for  the  self-consistent  peak 
field  amplitude  was  then  determined  by  involving  energy  balance:  i.e.,  the 
kinetic  energy  lost  by  electrons  transversing  the  cavity  is  equal  to  the 
microwave  power  emitted  from  the  cavity.  The  trajectories  of  a  set  of  test 
electrons  initially  distributed  randomly  with  respect  to  their  phase  angle 
in  velocity  space  were  integrated  numerically  across  the  cavity  in  the 
presence  of  the  electromagnetic  field  to  determine  the  average  fraction  of 
kinetic  energy  lost  to  the  field  for  a  given  peak  electric  field 
amplitude.  A  nonlinear  efficiency  function  was  constructed  which  gave  the 
fraction  of  beam  energy  converted  to  microwaves  as  a  function  of  peak 

electric  field  IeJ,  magnetic  field  B,  beam  voltage  VQ,  ratio  of  perpendiclar 
to  parallel  velocity^)  and  both  cavity  and  beam  geometric  factors.  The 
microwave  power  generated  by  the  beam  P^  was  given  by  the  following 


expression: 


THRESHOLD  CURRENT  VERSUS  MAGNETIC  FIELD 
FOR  DIFFERENT  INTERACTION  LENGTHS 


TABLE  I 


PARAMETERS  FOR  CARM  OSCILLATOR  DESIGN 


Beam  Power - - - 100  MW 

Beam  Voltage - 511  kV 

Beam  -  .  8 

Microwave  Power - ^30  MW 

Cavity  Mode - - -  TEqi 

Operating  Frequency  -  44  GHz 

Cavity  Cutoff  Frequency  -  10  GHz 

Cavity  Interaction  Length - 40  wavelengths 

Cavity  - - ^650 

Circuit  Magnetic  Field  - —  10300  gauss 


The  power  radiating  from  the  cavity  as  a  result  of  the  external  diffraction 


Q  was  given  by: 


«  (stored  energy) 


PR- 


kloss  | E  |  2 
■  o 


where  kloss  is  given  by: 


J^(k  )(1  -  a2/k2  )a2u_  zout  2 

kloss  -  (-S-S - - - 22 — S3.]  /  |E(z)|  dz 


zin 


Since  in  steady  state  Pr  «  Pg,  the  following  non-linear  equation  for  |Eo| 
was  obtained: 

n(|Eo|2,B,VQ,a,OVoIo  -  kloss  JeJ* 


In  practice,  the  current  IQ  was  determined  as  a  function  of  |Eep. 

Figure  3  shows  efficiency  versus  distance  for  an  optimal  electric  and 
magnetic  field  setting  indicating  that  over  35%  efficiency  could  be 
achieved.  The  higher  electric  field  was  optimal  for  a  shorter  length 
resonator.  Anticipating  that  non-ideal  beams  with  velocity  spread  would 
interact  more  efficiently  in  shorter  resonators,  a  detailed  efficiency 
analysis  was  performed  for  a  cavity  35  wavelengths  long.  Figure  4  shows 
power  out  versus  field  amplitude  squared  assuming  200  amps  of  current  for 
two  different  magnetic  fields. 

The  intersection  with  the  diffraction  Q  line  indicates  the  actual 
operating  points.  One  observes  that  the  magnetic  field  showing  higher 
efficiency  (dotted  curve)  has  two  solutions.  The  lower  power  solution  is  not 
physical  in  that  the  system  is  unstable  at  ths  point.  The  stable  high  power 
solution  would  appear  inaccessible  from  noise  (hard  excitation)  in  that  as 
E  -*■  o,  the  cavity  losses  are  greater  than  the  power  generated  by  the  beam. 
Accessibility  is  achieved  if  Q  is  increased  sufficiently  to  allow  the  Q 
line  to  have  a  smaller  slope  at  the  origin  than  the  power  curve.  For  this 


VOLTAGE  =  511  kV 


FIGURE  3  EFFICIENCY  VERSUS  DISTANCE  Z  FOR  OFHMAL  AND  NON-OPTIMAL 
ELECTRIC  FIELD  AMPLITUDES 
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increased  Q  value,  there  is  a  single  operating  point  intersection  but  at  a 
reduced  power  level.  In  practice,  the  finite  rise  time  characteristic  of 
the  beam  could  allow  accessibility  to  the  higher  power  states  without 
changing  the  Q  factor.  This  requires  further  investigation..  Figure  5 
shows  efficiency  versus  current  for  different  magnetic  fields:  The  dashed 
lines  are  the  non-physical  solutions.  An  efficiency  of  34%  is  indicated  at 
200  amps  resulting  in  a  predicted  output  power  of  34  MW.  The  optimum 
output  power  of  40  MW  is  illustrated  in  Figure  6,  where  a  3.0%  magnetic 
taper  was  applied  and  the  interaction  length  was  increased  to  40 
wavelengths . 

Section  III.  Preliminary  Designs  for  the  Magnetron  Injectin  Gun  and  Bragg 
Reflector 

A  preliminary  design  for  the  required  MIG  is  shown  in  Figure  7.  The 
initial  gun  was  designed  for  an  of  1.5  at  a  final  magnetic  field  of  5600 

gauss.  It  is  required  to  operate  the  gun  in  a  final  field  of  10300  gauss 
with  ano(  of  .8  and  at  the  same  guiding  center.  This  requires  that  the 
cathode  radius  be  changed  from  2.28  cm  to  3.  cm  and  the  cathode  magnetic 
field  be  increased  to  845  gauss.  Further,  refinements  of  the  electrode 
structure  would  be  required  to  ensure  beam  velocity  spreads  of <2%.  Table  II 
shows  the  design  parameters  of  the  modified  gun. 

Finally,  the  design  of  the  Bragg  reflector  is  required  which  gives  a  Q 
of  600.  In  a  single  mode  analysis  one  solves  the  following  differential 


equation  for-  the  complex  field  profile  E  and  complex  co : 

.?r/  x  oj2  _  u2  (z) 

d*E(z)  ,  r  cE  i  . 

~dP~~  I - Z2 - J  *  0 

where: 


u^Cz)  -  «*(z)  [l  -  m^CzJ/luan]  =  effective  cutoff  frequency 
2  2  2  2 

«d  (a)  -  c  k  /a  (z)  =  local  cutoff  frequency 
c  tnn 


FIGURE  5  EFFICIENCY  VERSUS  CURRENT  FOR  DIFFERENT  MAGNETIC  FIELDS 
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FIGURE  7  PRESENT  500  kV  200  A  inn  MT.G  GUN  TO  BE  MODIFIED  FOR  CARM 


TABLE  II 


GUN  PARAMETERS  FOR  CARM 


Beam  Voltage  - 

Anode  Voltage  - 

Beam  Current  - 

Final  Guiding  Center  - 

Final  Larmor  Radius  - 

Final  B^  - 

Final  B  ^ - 

Radius  Cathode  - 

Radius  Anode  - 

Axial  B-Field - 

Cathode  B  Field - 

Beam  Current/Space  Charge  Current  - 

Field  Gradients  at  Cathode - 


500 

kV 

150 

kV 

200 

Amps 

.877 

cm 

.1775 

cm 

.541 

cm 

.676 

cm 

3.0 

cm 

5.57 

cm 

10300 

gauss 

845 

gauss 

30 

% 

80 

kV/cx 

The  diffraction  Q  of  the  cavity  is  then  given  by: 

Q= 

An  approximate  expression  for  the  frequency  and  Q  is  given  by: 

U(L=  \l  COj+i  JW /lY1  C*'!*) 

CO g  ^ftL/c 
(/- 

where  R|  and  Rj  are  the  reflection  co-efficient  of  the  Bragg  Reflectors. 
Reflection  is  obtained  by  coherent  scattering  of  a  forward  wave  into  a 
backward  wave  which  occurs  when  the  ripple  wave  length  is  1/2  the  guide 
wavelength. 

The  mode  structure  of  Bragg  resonators  (cutoff  on  the  left  side)  with  Q 
factors  of  566  and  832  are  illustrated  in  Figures  8  and  9.  The  real  and 
imaginary  parts  (dotted  line)  of  the  field  profile  is  indicated. 
Physically,  the  imaginary  part  is  the  field  profile  90°  in  phase  from  the 
real  part. 

A  more  comprehensive  technique  for  evaluating  Bragg  resonators  is  in 
progress  using  scattering  matrix  theory  which  includes  all  TE  and  TM  modes 
capable  of  coupling  with  the  rippled  wall.  The  Bragg  reflectors  are 
characterized  by  scattering  matrices  [S^]  and  [S^]  which  describe  the 
reflected  amplitudes  of  a  set  of  modes  incident  on  the  reflector.  The 
eigenfrequency  and  Q  are  determined  from  the  following  determinant: 

Det  ||l  -  [Si]  [S2]||  =  0 

The  requirements  on  the  Bragg  reflector  necessary  to  maintain  mode  purity 
can  be  determined. 

Section  IV.  Future  Work 

The  previous  work  must  be  thought  of  as  a  preliminary  analysis.  Test 
particle  simulation  must  be  done  in  the  actual  resonator  fields  indicated 


I  =  566 
=  43.93  GHz 


Length  Straight  Section 
Length  Bragg  Reflector 
Ripple  Amplitude  aj/a0 


FIGURE  8  BRAGG  RESONATOR  FIELD  PROFILE 
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Q  =  832 
f  =  43.9  GHz 


Length  Straight  Section  =  70  Ag 

Length  Bragg  Reflector  =  85 

Ripple  amplitude  al/ao  = 


FIGURE  9  BRAGG  RESONATOR  FIELD  PROFILE 
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in  Figure  8  and  9.  Velocity  spreads  must  be  given  to  the  initial  particle 
distributions  to  evaluate  the  effect  of  beam  emittance.  Threshold  current 
analysis  must  be  done  for  competing  modes  especially  the  "gyrotron  modes" 
near  cutoff.  Finally,  a  fully  self-consistent  EM  particle  simulation 
including  space  charge  effects  is  required  to  verify  a  final  design. 
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CARM  OSCILLATOR  STUDY 


A  study  was  undertaken  to  design  and  analyze  a  CARM  oscillator  which 
could  operate  at  150  GHz  with  power  levels  in  the  mesawa‘t  range.  The 
vol+aCTe  was  chosen  to  be  400  kv  in  order  to  be  compatable  with  possible 
cv  applications  as  well  as  minimizing  power  supply  requirements.  The 
mode  of  operation  chosen  was  the  circularly  polarized  TE-^  implying  a 
Pierce/Wissler  conf iguration  for  beam  formation,  Beam  quality  require¬ 
ments  (velocity  spread)  limited  sun  pervience  to  Z  .2  or  beam  current 
£  50  amperes.  In  order  to  be  near  autoresonance,  the  phase  velocity 
was  chosen  to  be  2^  above  the  speed  of  lisht.  The  operation  frequency 
of  150  GHz  was  thus  5  times  above  cutoff  (29  GHz).  The  voltase 
restriction  of  400  kv  required  that  the  relativistic  cyclotron  frequency 
be  two  times  the  cutoff  frequency  (i.e.,  magnetic  field  *—  37000  gauss). 

An  interaction  diagram  is  shown  in  fisure  1. 

The  upper  intersection  point  is  the  CARM  while  the  lower  intersection 
point  represents  a  possible  backward  wave  oscillation,  at  ^3.0  GHz  ('“'4255 
a^ove  cutoff).  In  order  to  create  a  standing  wave  in  a  guide  far  above 
cutoff,  one  uses  a  Brass  reflector.  The  resonator  selects  out  the  axial 
mode  whose  wavelensth  is  double  the  Brans  ripple  length.  The  cavity  Q 
is  in  seneral  much  lower  for  all  other  resonant  wavelengths  (frequencies) 
includins  the  mode  near  cutoff  (J?*l).  This  is  illustrated  in  fisures 
2a  and  2b.  A  cavity  was  chosen  with  a  6  cm  interaction  ler.sth  (25 
wavelensths)  and  a  Brass  reflector  consistins  of  20  ripples  each  having 
an  amplitude  of  .014  cm.  The  operating  mode  at  150  GHz  then  had  an 
“^fecrive  S  number  of  50  and  an.  external  Q  of  B^3  •  The  J?=  1  mode  near 
cutoff  had  a  Q  of  214.  The  main  problem  with  this  configuration  is 
‘hat  there  could  be  mode  camoetition  with  the  backward  wave  oscillation 
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which  does  not  require  external  reflection  to  exist  (very  low  Q). 

Three  different  computer  models  were  developed  to  analyze  this 
oscillator  confisurat ion.  The  first  model  looks  at  the  CARM  interaction 
onlv.  It  assumes  that  onlv  a  single  mode  at  a  sinsle  frequency  exists, 
namely  the  resonant  mode  created  by  the  Brasg  reflector..  It  assumes 
that  th°  beam  does  not  modify  the  empty  cavity  field  profile.  Electrons 
are  integrated  borh  analytically  and  numerically  through  these  fields 
to  obtain  the  start  oscillation  currents  and  saturated  efficiencies 
respectively.  The  second  model  was  developed  to  look  at  interaction 
with  the  backward  wave  oscillation  only.  This  self-cons intant  simula¬ 
tion  model  is  based  on  the  slow  time  scale  orbit  equations  coupled  to 
Maxwell's  equations.  An  efficient  searching  scheme  is  used  to  find  the 
unknown  initial  magnetic  field  and  oscillation  frequency  which  allows 
+he  boundary  conditions  to  be  satisfied.  Finally  a  third  model  has 
been  developed  which  is  a  fully  time  dependent  particle  mode  which  can 
model  self-consistently  both  the  CARM  interaciion  and  backward  wave 
interaction  simultaneously,  ^he  code  couples  the  exact  time  dependent 

orbit  equations  (represented  bv  X  (t),V  (t)  )  with  the  followins  reduced 
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s®t  of  Maxwell's  equations  for  interaction  with  TE  modese 
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Figure  Sa  Field  profile  versus  distance  for  backward  wave 

interaction.  Frequency  =42 . 84GHz  .  Efficiency  = 4.45* 
Length  =6.3  cm.  Ream  currents  2  5  amps. 
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Figure  5b  Efficiency  versus  interaction  length  for  backward  wave 
oscillation.  Beam  currents  25  amps. 


iMv. 


HYBRID  CYCLOTRON  INSTABILITY  IN  A  FREE 
ELECTRON  LASER  WITH  AN  AXLAL 
MAGNETIC  FIELD 


T.  H.  Kho  and  A.  T.  Lin 

Physics  Department 
University  of  California 
Los  Angeles,  California  90024 

Received  April  20,  1987 


Th«  ixldl  guide  magnetic  field  In  a  free  electron  laser 
<FEL)  is  found  to  give  rise  to  a  "hybrid"  FEL-  cyclotron 
instability.  This  instability  occurs  at  a  modified  second 
harmonic  FEL  frequency  uhere  the  axial  ponderomotive  force 
Is  resonant  with  a  natural  frequency  of  electron  beam  axial 
motion  In  the  combined  vlggler  and  guide  magnetic  field 
High  efficiency.  In  excess  of  30  percent,  and  strong  growth 
rate  are  observed. 


I .  Introduction 

In  many  free  electron  laser  experiments,  an  axial 
magnetic  field  Is  added  to  guide  the  electron  beam  into  the 
interaction  region.1  This  axial  field  represents  a  new 
degree  of  freedom  of  the  system.  Previous  investigations 
on  the  impact  of  the  axial  field  on  wave-particle  inter¬ 
actions  have  concentrated  on  the  Instability  at  the  free 
electron  laser  resonance,  and  to  a  lesser  extent,  that  at 
the  cyclotron  maser  resonance. 2- 3. A 

In  this  paper,  we  present  results  from  computer  simu¬ 
lations  of  a  "hybrid"  FEL-cyclotron  instability  that  is 
found  to  occur  at  a  modified  second  harmonic  FEL  resonance. 
A  l-2/2  dimensional,  relativistic,  electromagnetic  particle 

S 

code  Is  used  to  simulate  the  Interaction  of  a  tenuous 
electron  beam  with  the  radiation  field  In  a  static  magnetic 
field  of  the  form, 
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^  Kho  and  Lin 

2  =  Bcz  ♦  By<xCosk„z  +  ySinkyZ)  <1) 

where  Bc  and  By  are  the  guide  and  wiggler  field  strengths, 
respectively,  and  \y=2n/kw  is  the  wiggler  period.  This 
wiggler  field  is  an  approximation  of  the  field  near  the 
axis  jOf  the  bifilar  helix  configuration  used  in  experi¬ 
ments,  and  is  accurate  for  kwr<<l,  where  r  Is  the  radius 
of  beam  excursion  about  the  axis.  The  equilibrium  steady 
state  orbits  for  this  magnetic  field  configuration  is  given 

by.6  ...  . 


where 


p  =  p0z  ♦  py(xCoskwz  +  ySinkyZ ) 


Pw~^wPo^^c”^wVo^  *  fly—eBy/ Ypmpc , 


Oc=eBc/Y0mec, 


and  y0  is  the  initial  Lorentz  factor  of  the  beam.  The 
usual  FEL  resonance  arises  from  the  transverse  oscillation 
which  induces  radiation  at  the  Doppler  shifted  wiggler 
frequency  u=<k+k  >v0.  To  understand  the  hybrid  instability, 
we  note  that  for  small  perturbations  about  the  steady  state 
in  the  combined  wiggler  and  axial  field,  the  electron  orbit 
is  stable  (in  the  absence  of  an  electromagnetic  field),  but 
now  contains  additional  modes  of  oscillation,  approximately 
given  by, 

px  =  PyCoskyZ  +  piCos(Dct+^)  -  qCos<2kyZ  -  Qct  +  ♦) 

Py  =  PySinkwz  *  PiSin(Qct+i)  -  qSfn(2kwz  -  0ct  +  ♦>  <3> 

Pz  =  p0  -  ByPiCosfkyZ  -  Oct  ♦  #> 


3w  =  Pw  /Pc 


1  02 

q  =  o  3W  Pi 


where  p^  and  i  are  constants  of  motion  and  Pj<<Py  has  been 
assumed.  There  is  now  a  Larmour  component  to  the  electron 
motion  and  the  axial  velocity  is  no  longer  constant,  but 
oscillates  at  the  cyclotron  and  wiggler  beat  frequency. 

The  third  term  in  eq.<3)  for  the  transverse  velocity 
arises  from  the  coupling  of  the  axial  oscillation  with  By, 
and  gives  a  modified  second  harmonic  FEL  resonance,  u  - 
<k+2kv)v0-  0C.  Electromagnetic  waves  with  frequencies  in 
the  neighbourhood  of  this  resonance  are  found  to  be 
unstable.  This  instability  arises  for  the  following  reason. 
An  EM  wave  with  frequency,  w=(k*2kv.)vp-  £lc .  will  produce  ar. 

e 

axial  ponderomot  1  ve  force  (-  -  vw  *  E^l,  which  is  resonant 
with  the  natural  frequency  of  beam  axial  oscillation. 
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Aa=lkuv0-  Qcl,  <eq.3>.  For  small  but  finite  Aa,  this 
resonant  ponderomotive  force  can  greatly  amplify  the  axial 
oscillation  In  the  electron  motion,  which,  through  Its 
cross  product  with  By,  drives  a  transverse  current  that 
reinforces  che  radiation  at  gj=<  k+2kw  >  Vp-Op  .  This  resonant 
amplification  of  electron  axial  oscillation  has  been 
observed  in  our  simulations. 


I I .  Simulation  Results 

In  order  to  make  the  effect  of  this  instability  more 
transparent,  simulations  are  performed  keeping  only  one 
k-mode  of  the  radiation  field  at  a  time.  This  is  equiva¬ 
lent  to  simulating  an  amplifier  system.  A  high  efficiency 
regime  appears  to  exist  for  EM  waves  that  are  simultaneous¬ 
ly  resonant  with  both  longitudinal  and  transverse  second 
harmonic  electron  oscillations.  This  occurs  when  <2kwv0- 
Qo^Qq-  kwv0)  or*  Qo/kyV0»1.5.  Results  for  this  case  are 
shown  below  for  the  following  parameters: 

yo=3.0  flw/kuvo=0.12  flc/kyV0=1.44 

vo/c=0 . 91  Vy/c=0.25 

Figure  1  shows  the  dispersion  diagram  for  EM  waves 
propagating  along  the  axial  magnetic  field  in  a  plasma  with 
drift  velocity  v0  Note  that  the  right  polarized  slow 
cyclotron  wave  is  resonant  with  the  uj=(k+2kw>v0-ftc  line  at 
k=6  and  with  cte  FEL  mode  w=(k+ky)v0  near  k=10  (k=k/ky). 

Not  shown  is  the  w=(k-  kw)v0+  Oc  line  which  is  parallel  to, 
and  Just  slightly  below  w=< k+2k„>v0-  Q,- ;  it  intersects  the 
cyclotron  wave  at  k=S.  The  right  polarized  fast  EM  wave  is 
nearly  coincident  with  the  Doppler-shifted  cyclotron 
resonance,  o^kVp-Klc,  for  all  values  of  k  shown. 

Exponential  growth  in  the  EM  wave  occurs  for  the  whole 
range  of  k  in  figure  i.  Maximum  growth  rate  occurs  at 
k= 6,  w=5.82  where  the  cyclotron  wave  and  u>=< k+2kw> v0-nc 
intersect  (flg.i),  with  a  value  of  0.076<dp  <w=w/ujp, 

Up= Annexing > ,  while  maximum  efficiency  occurs  at  ic=5. 

Figure  lb  shows  the  evolution  of  the  radiation  field  at  k=S 
and  k-10.  Figure  2a  3hows  the  power  spectrum  at  saturation 
for  k=5,  8,  and  10.  The  power  spectrum  for  k=5  shows  the 
excitation  of  the  right  polarized  cyclotron  wave  by  the 
hvbrid  FEL-  cyclotron  resonance.  For  k=8,  two  instabilities 
are  present.  The  lower  frequency  peak  in  the  power 

spectrum  corresponds  to  the  right  polarized  slow  cyclotron  • 
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Figure  1.  Dispersion  diagram  for  G0=1.28,  yo=3.0  and  3o=0.91. 

Lines  shown:  (k+kw)v0  (1);  u>=  < k+2ku> v0-  00  (2); 

right  polarized  fast  EM  <3>,  and  slow  cyclotron 
E(1  (unlabelled)  waves. 


wave  excited  by  the  hybrid  resonance.  The  higher  frequency 
peak  is  the  FEL  mode  w=<k+kv)v0.  For  k-10,  only  one  mode  is 
evident.  This  is  at  the  FEL  resonance,  and  the  growth 
rate,  at  0.037up,  is  about  half  that  of  the  hyorld  instabi¬ 
lity  (fig.  2b).  Conspicuously  absent  from  these  simula¬ 
tions  is  anv  sign  of  an  Instability  at  the  Doppler-shifted 
cyclotron  resonance,  u>=kv0+fiCi  even  though  the  right 
polarized  fast  EM  wave  frequency  is  very  close  to  this 
resonance  for  the  values  of  k  here. 


In  an  FEL  instability,  the  electron  beam  Is  bunched 
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axially  by  the  ponderomottve  force  at  a  periodicity  of  k+kw 
land  its  harmonics).  This  occurs  here,  but  in  addition.  3 
well  defined  change  in  the  pitch  of  the  electron  beam  is 
observed  in  the  evolution  of  the  second  harmonic 
instability  at  k-S  land  also  at  k--6) .  Ftgnre  3  displays 
the  azimuthal  location  of  electrons  about  the  axis, 
parameterized  by  B^tan-1 1  Pvi/Pzi  >  ■  over  one  wiggler 
period,  at  the  initial  time  (a),  at  Upt=l40  (just  before 
saturation)  for  k=S  (b),  and  at  Wpt=2S0  for  k-10  <c>.  The 
electron  beam,  which  initially  is  a  helix  with  the  same 
pitch  as  the  magnetic  wlggler,  evolves,  in  the  case  of  k-S, 
into  a  helix  with  a  pitch  of  k+2kw.  This  new  pitch  results 
from  the  3trong  transverse  Lorentz  force  on  the  beam, 
provided  by  Bw  and  the  large  axial  oscillation  In  velocity 
driven  bv  the  resonant  ponderomotive  force.  At  k-10.  awav 
from  the  hybrid  FEL-cyclotron  resonance,  the  electron  beam 
preserves  its  initial  helicity  throughout  the  unstable  wave 
growth.  Electron  bunching  at  k*kw  can  be  seen  for  both 
k=S  and  k=10. 

Maximum  efficiency  is  obtained  at  k-S,  where  3S  percent 
of  beam  kinetic  energy  Is  converted  into  radiation  energy. 
This  efficiency  is  high,  although  it  is  only  about  half  of 
what  might  have  been  expected  from  the  reduction  in  beam 
axial  velocity  at  saturation,  vs/vo»0.88  < vs-vz(Upt -ISO) > . 

If  all  the  loss  in  axial  kinetic  energy  goes  into 
radiation,  the  efficiency  can  be  estimated  from  n=(y0- 
rs)/(r0~l  )s,0- 7  (ys-1/^(  1~3s>  .  3s-vg/c>.  The  reason  for 
this  discrepancy  is  that  a  large  fraction  of  the  axial 
kinetic  energy  has  been  converted  into  azimuthal  energy. 

The  mean  transverse  velocity  has  increased  by  about  a 
factor  of  two  at  saturation  over  its  initial  value. 

As  noted  above,  at  k- 5,  the  cyclotron  wave  is  resonant 
with  the  axial  oscillation  of  the  beam  leq.  3>.  This 
appears  to  be  relevant  to  the  achivement  of  high  efficiency 
here.  Figure  4  shows  plots  of  kinetic  energy  <y-l>  versus 
axial  velocity  of  electrons  at  saturation,  upt-150,  and 
after  saturation  at  uipt=270,  for  k=5.  The  electrons  are 
initially  concentrated  at  the  point  <r-l>=2  and 
3z^vz/c=0.91 .  The  resonance  line  for  u>=<k-kw>vz+0/r  is 
shown  in  figure  4a,  for  k-- 5,  w-4.88  (the  observed  radiation 
frequency),  where  fi=eB0/mec.  Electrons  on  this  line  have 
axial  oscillations  that  are  resonant  with  the  EM  wave.  At 
saturation,  the  bulk  of  the  electrons  have  been  shifted 
down  into  the  !y~1)  vs.  vz  phase  space  encompassing  this 

7 

line,  suggesting  some  degree  of  autoresonance.  For  the 
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FEL  resonance  at  k=10,  the  efficiency  Is  only  five  percent 
and  the  <y-l)  vs.  vz  plot  (figure  4c)  does  not  display  the 
above  feature. 

After  the  hybrid  Instability  saturated,  some  electrons 
continue  to  be  accelerated  by  the  EM  wave  to  high  energies 
(figure  4b>.  There  are  two  distlnguislble  groups  of 
electrons;  one  supporting  the  hybrid  wave,  the  other  being 
accelerated  by  it.  The  accelerated  electrons  appear  to  be 
confined  to  a  fin  shaped  region  In  (y-1)  vs.  vz  space,  and 
bound  at  the  upper  vz  limit  by  the  cyclotron  autoresonance 
line,  shown  in  figure  4b  (phase  velocity  of  EM 

wave  at  k-S  is 'w/kc =0.998) . 


III.  Conclusions 

The  introduction  of  an  axial  magnetic  field  into  a 
helical  wiggler  field  introduces  the  possibility  of  a 
natural  low  frequency  mode  of  axial  oscillation  for  the 
electron  beam.  When  driven  by  a  resonant  ponderomotlve 
force  associated  with  an  EM  wave  at  the  hybrid  FEL- 
cyclotron  resonance,  w=(k*2kw)v0-  Clc,  this  axial 
oscillation  can  reach  large  amplitude  so  that  when 
translated  into  transverse  oscillation  by  Bw,  it  serves  to 
amplify  the  radiation.  It  is  shown  that  some  degree  of 
autoresonance  operate  when  the  EM  wave  Is  also  approximate¬ 
ly  resonant  with  the  axial  oscillation,  enabling  high 
efficiency,  in  this  example,  In  excess  of  30  percent,  to  be 
achieved.  In  a  one-dimensional  plasma,  this  (empirical) 
condition  <nc/kvv0-l . S)  restricts  the  hybrid  resonance  to 
the  slow  cyclotron  wave.  However,  this  restriction  is  not 
necessary  in  a  waveguide  configuration,  where,  for 
sufficiently  large  nc  the  hybrid  resonance  can  occur  on  a 
fast  EM  wave.  The  efficiency  here  compares  very  favourably 
with  a  cyclotron  maser  operated  near  autoresonance.  For  a 
Lorentz  factor  of  Y0=3,  the  expected  efficiency  of  a 
cyclotron  maser  is  about  20  percent,  and  for  a  relativistic 

Q 

gvrotron  .  is  about  IS  percent.  However,  to  obtain  optimum 
efficiency,  the  electron  beam  has  to  be  injected  with  high 
uniformity  In  azimuthal  velocity,  which  is  technically 
difficult  to  achieve.  Furthermore,  the  gain  near 
autoresonance  (w/ksc)  In  a  cyclotron  maser  is  very  small. 

In  contrast,  for  the  hybrid  instability,  no  prior  tailoring 
of  beam  transverse  velocity  is  needed,  and  the  gain  is 
relatively  high.  These  advantages  makes  the  hybrid 
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Instability  a  potentially  attractive  alternative  to  the 
cyclotron  aaser  as  a  high  efficiency  microwave  source. 
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